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Abstract:  The  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory  is  deveioping  a  mechanistic 
pavement  design  procedure  for  use  in  seasonai  frost 
areas.  The  procedure  was  used  to  predict  pavement 
performance  of  some  test  sections  under  construction 
at  the  Mn/ROAD  faciiity.  Simulations  were  conducted  in 
three  phases,  investigating  the  effects  on  predictions  of 
water  table  position,  subgrade  characteristics,  asphait 
model,  and  freeze  season  characteristics.  The  proce¬ 
dure  predicted  significantiy  different  performance  by 
the  different  test  sections  and  highiy  variable  results 
depending  on  the  performance  model  applied.  The 
simulated  performance  of  the  tests  sections  also  was 


greatly  affected  by  the  subgrade  conditions,  e.g.,  den¬ 
sity,  soil  moisture,  and  water  table  depth.  In  general, 
predictions  for  the  full-depth  asphalt  sections  indicate 
that  they  will  not  fail  due  to  cracking,  but  two  of  the 
three  criteria  for  subgrade  rutting  indicate  failure  before 
the  five-  or  1 0-year  design  life  of  the  sections.  Conven¬ 
tional  sections  are  predicted  not  to  fail  due  to  subgrade 
rutting;  however,  sections  including  the  more  frost- 
susceptible  bases  in  their  design  are  predicted  to  fail 
due  to  asphalt  cracking  relatively  early  in  their  design 
life,  and  sections  with  non-frost-susceptible  bases  are 
predicted  to  fail  towards  the  end  of  the  design  life. 
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EXECUTIVE  SUMMARY 


The  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  (CRREL)  is 
developing  a  mechanistic  pavement  design  procedure  for  use  in  seasonal  frost  areas. 
It  consists  of  four  computer  programs  that  compute  soil  and  pavement  moisture  and 
temperature  conditions  (FROST),  resilient  modulus  and  Poisson’s  ratio  (TRANS¬ 
FORM),  stresses  and  strains  in  the  pavement  system  (NELAPAV),  and  cumulative 
damage  (CUMDAM).  Damage  predictions  are  based  on  several  equations  that 
employ  horizontal  strain  at  the  base  of  an  asphalt  layer,  vertical  strain  at  the  top  of 
the  subgrade,  or  horizontal  stress  at  the  base  of  a  concrete  layer.  The  procedure  was 
used  to  predict  pavement  performance  of  some  of  the  Mn/ROAD  test  sections. 
Laboratory  tests  on  the  Mn/ROAD  materials  (Bigl  and  Berg  1996a,  Berg  et  al. 
1996)  provided  the  input  parameters  necessary  for  the  modeling  effort.  Simulations 
were  conducted  in  three  phases,  investigating  the  effects  on  predictions  of  water 
table  position,  subgrade  characteristics,  asphalt  model,  and  freeze  season  character¬ 
istics. 

Phase  1,  conducted  in  the  spring  of  1991,  included  an  initial  simulation  series 
that  modeled  temperatures  from  a  year  close  to  the  mean  freezing  index.  These 
boundary  conditions  were  applied  to  eight  flexible  and  three  rigid  sections.  Phase  2, 
an  effort  in  the  summer  of  1992,  had  two  primary  objectives.  Phase  2A  included 
three  series  modeling  the  eight  flexible  sections  with  the  mean  freeze  season  and 
changing  the  method  employed  to  calculate  the  asphalt  and  subgrade  moduli.  Phase 
2B  investigated  the  variability  in  predictions  when  temperatures  from  freeze  sea¬ 
sons  with  maximum  and  minimum  freezing  indices  are  applied  to  a  single  flexible 
section.  Phase  3,  an  effort  in  the  summer  of  1993,  expanded  the  investigation  of  the 
effects  of  freeze  season  characteristics.  This  series  modeled  21  different  freeze 
seasons  applied  to  one  full-depth  and  one  conventional  flexible  section. 

The  procedure  predicted  significantly  different  performance  by  the  different  test 
sections  and  highly  variable  results  depending  on  the  performance  model  applied. 
The  simulated  performance  of  the  test  sections  was  also  greatly  affected  by  the 
subgrade  conditions,  e.g.,  density,  soil  moisture  and  water  table  depth.  In  general, 
predictions  for  the  full-depth  asphalt  sections  indicate  that  they  will  not  fail  due  to 
cracking,  but  two  of  the  three  criteria  for  subgrade  rutting  indicate  failure  prior  to 
the  5-  or  10-year  design  life  of  the  sections.  Conventional  sections  are  predicted  not 
to  fail  due  to  subgrade  rutting,  but  sections  including  the  more  frost-susceptible 
bases  in  their  design  are  predicted  to  fail  due  to  asphalt  cracking  relatively  early  in 
their  design  life,  and  sections  with  non-frost-susceptible  bases  are  predicted  to  fail 
towards  the  end  of  their  design  life. 
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Modeling  of  Mn/ROAD  Test  Sections 
with  the  CRREL  Mechanistic 
Pavement  Design  Procedure 

SUSAN  R.  BIGL  AND  RICHARD  L.  BERG 


INTRODUCTION 

This  is  one  of  four  reports  related  to  the  full- 
scale  pavement  test  facility,  Mn/ROAD,  con¬ 
structed  by  the  Minnesota  Department  of  Trans¬ 
portation  (Mn/DOT)  adjacent  to  Interstate  94  in 
Otsego,  Wright  County.  One  report  describes  the 
results  of  laboratory  tests  conducted  to  determine 
the  physical  and  freeze-thaw  related  characteris¬ 
tics  of  the  on-site  subgrade  and  two  of  the  pave¬ 
ment  system  materials  used  in  the  test  sections 
(Bigl  and  Berg  1996a).  Another  report  discusses 
resilient  modulus  tests  conducted  on  the  above 
materials  (Berg  et  al.  1996).  This  report  provides 
an  initial  description  of  the  computer  programs 
used  in  a  mechanistic  pavement  design  procedure 
under  development  by  CRREL,  and  then  describes 
a  computer  modeling  effort  that  utilizes  the  re¬ 
sults  of  the  material  testing  program  in  the  proce¬ 
dure  to  predict  the  performance  of  some  Mn/ 
ROAD  test  sections.  The  final  report  of  the  series 
summarizes  information  in  the  first  three  reports 
(Bigl  and  Berg  1996b). 

The  mechanistic  pavement  design  procedure  is 
being  developed  for  use  in  seasonal  frost  areas 
and  considers  seasonal  variations  in  pavement 
strength  such  as: 

1 .  Large  increases  in  base  and  subgrade  strengths 
when  frozen. 

2.  Loss  of  base  and  subgrade  strengths  during 
the  spring  when  thawing  and  excess  water 
weakens  the  layers. 

3.  Loss  of  asphaltic  concrete  strength  during 
the  summer  months  when  asphalt  moduli 
decrease  with  higher  temperatures. 


Availability  of  data  from  the  Mn/ROAD  facil¬ 
ity  will  allow  us  to  verify  and  refine  various  as¬ 
pects  of  the  model. 

COMPUTER  MODELS 

The  mechanistic  design  procedure  consists  of 
the  four  programs  FROST,  TRANSFORM, 
NELAPAV,  and  CUMDAM  (Fig.  1).  FROST  pre¬ 
dicts  the  amount  of  frost  heave  and  thaw  settle¬ 
ment  of  the  pavement  structure  and  conditions 
throughout  the  depth  of  the  structure  (tempera¬ 
ture,  water  content,  pore  water  pressure,  ice  con¬ 
tent,  density)  at  a  given  time  increment.  TRANS¬ 
FORM  uses  the  output  from  FROST  as  input  and 
divides  the  pavement  structure  into  “layers”  based 
on  moisture  or  temperature  conditions.  Each  layer 
is  then  assigned  a  resilient  modulus,  Poisson’s 
ratio  and  density  value.  NELAPAV,  a  nonlinear 
layered  elastic  program,  calculates  stresses,  strains 
and  deflections  at  specified  locations  within  a  pave¬ 
ment  profile  when  a  load  is  applied  to  the  surface. 
CUMDAM  calculates  the  incremental  and  cumu¬ 
lative  damage  the  pavement  undergoes  using  sev¬ 
eral  available  damage  models.  The  behavior  is 
normally  modeled  for  a  one-year  period,  with  out¬ 
put  on  a  daily  basis.  By  assuming  that  the  results 
from  a  single  year  will  be  repeated  annually,  the 
number  of  applications  to  failure  is  estimated. 

FROST 

FROST  is  a  one-dimensional  coupled  heat  flow 
and  moisture  flow  model  that  computes  frost  heave 
and  thaw  settlement  of  a  pavement  or  soil  profile 


Output/Input 


with  time.  It  also  calculates  temperature,  moisture 
stress,  water  content,  ice  content,  and  density 
through  the  depth  of  the  profile  at  each  time  incre¬ 
ment.  FROST  was  originally  developed  by  Berg 
et  al.  (1980)  in  a  cooperative  study  funded  by  the 
Corps  of  Engineers,  the  Federal  Highway  Admin¬ 
istration,  and  the  Federal  Aviation  Administra¬ 
tion.  Additional  details  beyond  those  described 
here  are  given  in  Guymon  et  al.  (1993).  The  model 
assumes  one-dimensional  vertical  heat  and  mois¬ 
ture  flux  and  is  based  on  a  numerical  solution 
technique  termed  the  nodal  domain  integration 
method.  The  nodal  domain  integration  method  al¬ 
lows  use  of  the  same  computer  program  to  solve  a 
problem  by  either  the  finite  element  method,  the 
integrated  finite  difference  method,  or  any  other 
mass  lumping  numerical  method.  For  this  study, 
we  chose  to  use  the  integrated  finite  difference 
computational  method. 

Figure  2  shows  how  FROST  uses  nodes,  which 
are  exact  points,  to  divide  the  column  of  material 
into  horizontal  elements.  Material  properties  are 
assigned  to  each  element. 

The  program  was  developed  for  use  on  prob¬ 
lems  of  seasonal  freezing  and  thawing  of  non¬ 
plastic  soils,  and  is  based  on  the  following  pri¬ 


mary  assumptions,  with  additional  assumptions 
reported  in  Guymon  et  al.  (1993): 

1.  Darcy’s  law  applies  to  moisture  movement 
in  both  saturated  and  unsaturated  conditions. 

Element  Node 


Actual  Soli  Column  Mathematical 

Representation  Model  Column 


Figure  2.  Example  of  pavement  profile  divided 
into  finite  elements. 
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2.  The  porous  media  are  nondeformable  as  far 
as  moisture  flux  is  concerned;  i.e.,  consoli¬ 
dation  is  negligible. 

3.  All  processes  are  single  valued;  i.e.,  hyster¬ 
esis  is  not  present  in  relationships  such  as 
the  soil  water  characteristic  curve. 

4.  Water  flux  is  primarily  as  liquid;  i.e.,  vapor 
flux  is  negligible. 

The  governing  equation  used  in  FROST  to  de¬ 
scribe  soil  moisture  flow  is  derived  by  substitut¬ 
ing  the  extended  Darcy  moisture-flow  law  into 
the  one-dimensional  continuity  equation  for  an 
incompressible  fluid  flowing  through  porous 
media,  i.e.. 


dx 


^  I  Pi 
dt 


96i 


(1) 


where  = 
h  = 

X  = 
%  = 
Pi  = 
Pw  = 

0i  = 

t  = 


unsaturated  hydraulic  conductivity 
(permeability)  (cm/hr) 
total  hydraulic  head  (cm  water) 
depth  (cm) 

volumetric  unfrozen  water  content  (%) 
density  of  ice  (g/cm^) 
density  of  water  (g/cm^) 
volumetric  ice  content  (%) 
time  (hr). 


The  total  hydraulic  head  h  equals  the  sum  of  the 
pore  pressure  head  {hp  =  u/p^),  where  u  is  the  pore 
water  pressure,  plus  the  elevation  head  {h^  =  -x), 
where  x  is  measured  vertically  downward.  The  ice 
sink  term,  pj  30j  /p^  dt,  exists  only  in  freezing  or 
thawing  zones,  and  in  these  zones,  eq  1  is  coupled 
to  the  heat  transport  equation.  The  ice  sink  term 
assumes  that  0j  is  a  continuous  function  of  time. 

In  FROST,  the  soil  water  characteristics  are 
represented  using  a  relationship  in  the  form  of 
Gardner’s  (1958)  equation: 


0 


u  ~ 


(2) 


For  each  soil  to  be  modeled,  point  values  of  0u 
and  hp  are  determined  in  a  laboratory  moisture 
retention  test  (Ingersoll  1981).  Gardner’s  eq  2  is 
then  fit  to  the  data  using  a  least  squares  approach 
to  determine  the  best  fit  parameters  A,^  and  a.  Test 
results  for  the  Mn/DOT  materials  are  given  in  the 
first  report  of  this  series  (Bigl  and  Berg  1996a). 

Unsaturated  hydraulic  conductivity  is  also  ap¬ 
proximated  in  FROST  using  a  Gardner’s  equa¬ 
tion: 


(3) 


where 

hp 

P 


unsaturated  hydraulic  conductivity 
(cm/hr) 

saturated  hydraulic  conductivity 
(cm/hr) 

pore  pressure  head  (cm  of  water) 
Gardner’s  multiplier  for  hydraulic 
conductivity 

Gardner’s  exponent  for  hydraulic 
conductivity. 


Point  values  of  and  hp  for  each  soil  are  deter¬ 
mined  in  the  laboratory  by  an  unsaturated  hydrau¬ 
lic  conductivity  test  (Ingersoll  1981),  and,  again, 
Gardner’s  eq  3  is  fit  to  the  data  using  a  least 
squares  approach  to  determine  the  best  fit  param¬ 
eters  A^and  p.  See  Bigl  and  Berg  (1996a)  for  the 
test  results  on  the  Mn/ROAD  materials. 

Within  the  partially  frozen  zone,  FROST  re¬ 
duces  the  unsaturated  hydraulic  conductivity  us¬ 
ing  an  empirical  constant,  termed  the  E-factor, 
combined  with  the  ice  content  according  to  the 
following  equation: 

A:p  =  A:H(/!p)xio-'^®-,Eei>o  (4) 

where  ATp  =  the  adjusted  hydraulic  conductivity 
in  a  partially  frozen  element  (cm/hr) 
£*  =  an  empirical  constant,  dimensionless 
0j  =  volumetric  ice  content  (%). 


where  0^  =  volumetric  unfrozen  water  content  (%) 
0Q  =  soil  porosity  (%) 

/ip  =  pore  pressure  head  (cm  of  water) 

=  Gardner’s  multiplier  for  the  mois¬ 
ture  characteristics 

a  =  Gardner’s  exponent  for  the  moisture 
characteristics. 


For  this  study,  the  f'-factor  was  determined  within 
the  FROST  program  using  an  empirically  derived 
equation  based  on  the  saturated  hydraulic  conduc¬ 
tivity,  k^,  in  centimeters/hour: 

£  =  l(/:^_3f+6.  (5) 
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Frost  heave  is  estimated  from  the  total  amount 
of  ice  segregation  in  the  frozen  zone  by: 

0$  =  0i  “  (Qq  ~  ®n) 

where  0^  =  volumetric  segregated  ice  content  (%) 
00  =  porosity  (%) 

00  =  residual  unfrozen  water  content  (%). 

If  05  is  greater  than  0,  ice  segregation  has  occurred 
and  the  frost  heave  is  computed  by  multiplying  0s 
by  the  zone  thickness.  The  0^  parameter  estab¬ 
lishes  the  pore  water  stress  at  the  freezing  front 
for  the  solution  of  the  moisture  transport  equa¬ 
tion.  In  this  study,  0n  was  obtained  by  assuming  a 
moisture  tension  of  -800  cm  of  water  and  solving 
eq  2.  The  use  of  the  -800  cm  of  water  condition 
stems  from  that  being  the  highest  tension  mea¬ 
sured  in  various  field  studies.  Thaw  settlement 
from  ice  melting  is  the  reverse  process  of  that 
described  above  for  ice  segregation. 

To  conduct  the  calculations  described  above, 
FROST  requires  the  following  input  for  each  ma¬ 
terial:  1)  Gardner’s  coefficients  for  soil  moisture 
characteristics,  2)  Gardner’s  coefficients  for  hy¬ 
draulic  conductivity  characteristics,  3)  porosity  and 
density  of  the  soil,  4)  thermal  conductivity  and 
volumetric  heat  capacity  of  the  dry  soil,  and  5)  the 
£-factor. 

FROST  also  requires  the  following  input  for 
initial  and  boundary  conditions:  1)  element  lengths, 
2)  upper-  and  lower-boundary  pore  water  pres¬ 
sures,  3)  upper-  and  lower-boundary  temperatures, 
4)  initial  temperature,  pore  pressure  and  ice  con¬ 
tent  distributions  with  depth,  5)  surcharge  pres¬ 
sure,  6)  freezing  point  depression  and  7)  modifier 
of  the  upper  node  during  thaw. 

In  all  cases,  the  pavement  structure  was  simu¬ 
lated  as  a  column  with  its  upper  boundary  at  the 
pavement  surface  and  extending  down  to  400  cm 
(13.1  ft)  using  99  elements.  The  length  of  ele¬ 
ments  within  the  expected  zone  of  freezing 
(down  to  110  cm  or  3.6  ft)  was  about  2  cm  (0.8 
in.).  These  lengths  were  adjusted  for  individual 
cases  to  provide  nodes  positioned  exactly  at  the 
depths  where  the  interface  between  materials 
are  located  in  the  Mn/ROAD  test  sections.  The 
deeper  soil  was  modeled  with  element  lengths 
of  4  cm  between  1 10  and  230  cm  (3.6-7.6  ft),  10 
cm  between  230  and  340  cm  (7.6— 1 1.2  ft),  and 


20  cm  from  340  cm  (11.2  ft)  to  400  cm  (13.1  ft). 

The  upper  boundary  pore  water  pressure  was 
chosen  to  be  computer-generated,  as  follows. 
When  the  profile  is  completely  thawed  and  down¬ 
ward  vertical  drainage  occurs,  the  surface  pore 
water  boundary  condition  is  modeled  by 


which  means  that  the  velocity  flux  across  this 
boundary  is  zero.  The  upper-boundary  condition 
is  set  to  0  cm  of  water  when  the  upper-boundary 
temperature  is  above  0°C  and  frozen  regions  re¬ 
main  in  the  column.  When  the  surface  tempera¬ 
ture  is  below  0°C,  a  specified  constant  upper¬ 
boundary  pore  pressure  is  used.  To  be  consistent 
with  previous  studies,  a  value  of -300  cm  of  water 
was  used. 

The  lower-boundary  pore  pressure  condition  of 
FROST  is  set  by  specifying  discrete  pore  water 
pressures  that  relate  to  the  water  table  elevation  at 
times  when  these  conditions  occur.  At  intermedi¬ 
ate  times,  lower-boundary  pore  water  pressures 
are  linearly  interpolated.  For  all  the  Mn/DOT  cases, 
we  set  the  lower  boundary  pore  pressure  to  pro¬ 
duce  a  constant  water  table  depth  throughout  the 
simulation.  We  simulated  the  water  table  in  each 
test  section  at  the  depth  determined  by  field  mea¬ 
surements  to  be  representative  of  the  on-site  con¬ 
ditions.  Where  the  measured  water  table  varied 
through  a  test  section,  we  conducted  two  simula¬ 
tions  using  the  deepest  and  shallowest  values. 

Input  for  the  upper  boundary  temperature  con¬ 
dition  consists  of  a  set  of  specified  times  and  tem¬ 
peratures  that  are  implemented  as  step  changes. 
Values  were  input  in  24-hr  increments  using  the 
mean  daily  air  temperature.  Conditions  at  Buf¬ 
falo,  Minnesota,  were  simulated,  since  this  is  the 
nearest  station  to  the  Mn/ROAD  facility  (16  km; 
10  mi)  with  a  reasonably  long  record  of  meteoro¬ 
logical  data.  The  time  period  simulated  was  1 
October  1959  to  14  November  1960.  If  the  sever¬ 
ity  of  a  winter  is  judged  by  its  air  freezing  index, 
the  1959-1960  winter  is  very  near  the  average 
value  for  the  28-year  period  ending  in  1987.  The 
distribution  of  freezing  indices  at  Buffalo  during 
this  time  is  shown  in  Figure  3.  Starting  the  simu¬ 
lations  on  1  October  gave  30  days  of  computa¬ 
tions  before  the  first  freeze  event,  allowing  the 
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Figure  3.  Distribution  of  seasonal  freezing  index  with  time  at  Buffalo, 
Minnesota. 


model  to  stabilize  and  gradational  moisture  and 
temperature  conditions  to  develop. 

FROST  adjusts  the  air  temperature  values  to 
represent  the  soil  surface  temperatures  using  a 
procedure  similar  to  the  Corps  of  Engineers  n- 
factor  approach  for  seasonal  freezing  indices  (U.S. 
Army  1966).  An  n-factor  is  defined  as  the  ratio  of 
the  surface  index  to  the  air  index,  separately  cal¬ 
culated  for  the  full  freeze  and  thaw  season.  In  this 
study,  0.5  was  used  for  the  freezing  n-factor  (i.e., 
on  days  when  the  air  temperature  was  <  0°C),  and 
1 .8  was  used  for  the  thawing  n-factor. 

Bottom  boundary  temperature  conditions  consist 
of  a  set  of  times  and  temperatures  that  are  linearly 
interpolated  at  intermediate  times.  Three  points 
were  set  during  the  simulations  using  estimated 
values  of  a  reasonable  ground  temperature.  The  start 
and  end  temperatures  are  estimated  to  be  9.59°C 
(49.3°F),  which  is  3°C  (5°F)  higher  than  the  mean 
annual  air  temperature  6.83°C  (44.3°F).  The  third 
point  was  set  so  that  a  minimum  value  of  6.0°C 
(42.8°F)  would  occur  at  the  end  of  the  freeze  season. 

Initial  conditions  required  to  be  set  for  the 
FROST  program  are  temperature,  moisture  stress, 
and  ice  content.  The  initial  temperature  of  all  ele¬ 
ments  was  estimated  to  be  5°F  higher  than  the 
mean  annual  air  temperature.  For  cases  with  wa¬ 
ter  table  depths  less  than  3  m  (10  ft),  the  initial 
moisture  tension/pressure  was  set  at  a  gradient  to 
produce  a  water  table  (i.e.,  pressure  =  0)  at  the 
specified  depth,  increasing  positively  downwards 
at  1  cm  water/1  cm  depth.  For  simulations  with 


intermediate  water  table  depths  (3.6,  5.6  and  6.1 
m;  12, 18,  and  20  ft),  the  initial  tension  was  set  to 
-200  cm  water  (an  estimate  of  optimum  condi¬ 
tions)  from  the  surface  to  the  depth  beneath  which 
the  tension  from  the  theoretical  gradient  would 
increase  positively.  In  simulations  with  the  deep¬ 
est  in-situ  water  tables  (13.7  and  15.2  m  or  45  and 
50  ft),  we  set  a  minimum  bottom  tension  of  -300 
cm,  with  tensions  decreasing  from  the  bottom  to  2 
m  (6.5  ft),  and  a  constant  tension  of  -200  cm 
water  from  2  m  (6.5  ft)  to  the  surface.  Initial  volu¬ 
metric  ice  content  was  set  at  0.0%  for  all  ele¬ 
ments,  since  the  simulations  began  prior  to  the 
start  of  the  freeze  season. 

Surcharge  (overburden)  pressure  is  a  constant 
value  to  simulate  the  pressure  acting  on  the  top 
node  of  the  modeled  column.  In  a  case  where  only 
soil  beneath  a  pavement  surface  is  being  modeled, 
this  overburden  should  simulate  the  pressure  of 
the  pavement  on  the  soil.  For  the  Mn/ROAD  test 
section  simulations,  the  pavement  properties  were 
included  in  the  modeled  column,  so  the  overbur¬ 
den  pressure  was  set  at  zero. 

The  freezing  point  depression  is  a  constant  value 
that  represents  the  temperature  at  which  in-situ  soil 
water  freezes.  This  was  set  in  all  cases  to  be  0°C. 

A  modifier  is  included  that  adjusts  the  overbur¬ 
den  pressure  acting  on  the  upper  node  during  thaw 
periods.  This  value  can  be  set  between  0.0,  to 
represent  old,  cracked  pavements,  and  1 .0,  to  rep¬ 
resent  brand  new  pavements.  We  used  1.0  for  the 
Mn/ROAD  cases. 
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Standard  procedures  of  FROST  calculations 
were  selected  as  follows.  The  fully  implicit  method 
was  used  for  the  moisture  time  domain  solution, 
and  the  Crank-Nicolson  method  for  the  heat  trans¬ 
fer  time  domain  solution.  Simulations  were  mn 
with  a  time  step  of  0.2  hours,  which  is  the  time  at 
which  boundary  conditions  are  adjusted;  updates 
of  the  thermal  and  hydraulic  properties  were  set  to 
occur  once  per  hour. 

FROST  produces  what  can  be  a  very  large  out¬ 
put  file.  The  first  part  of  the  file  is  a  listing  of  the 
initial  profile  conditions,  including  all  the  input 
material  properties.  The  next  section  is  an  incre¬ 
mental  (usually  daily)  listing  of  the  conditions 
generated  for  each  node,  including  temperature, 
pore  water  pressure,  water  content,  ice  content, 
density,  and  porosity.  The  final  section  of  the  out¬ 
put  file  is  a  summary  of  the  predicted  frost  heave 
and  frost/thaw  penetration  data. 

TRANSFORM 

TRANSFORM  was  developed  at  CRREL  by 
Chamberlain  et  al.  (in  prep.)  and  was  modified 
extensively  by  Wendy  Allen  and  Gregor  Fellors, 
both  of  CRREL,  for  this  study.  The  TRANSFORM 
program  uses  FROST  daily  output  files  as  input 
and  produces  files  of  layered  pavement  systems 
where  each  layer  is  assigned  a  resilient  modulus, 
Poisson’s  ratio,  density  and  thickness.  Output  files 
from  TRANSFORM  are  in  the  format  to  be  used 
as  input  files  to  NELAPAV,  the  layered  elastic 
program  used  to  compute  stresses  and  strains  in 
the  pavement  system.  NELAPAV  requires  a  single 
file  containing  both  load  and  stmctural  data  for 
each  day  of  the  simulation. 

TRANSFORM  first  reads  from  a  separate  file  a 
series  of  material  type  identifiers  and  their  associ¬ 
ated  parameters  needed  for  calculating  the  moduli 
of  the  materials  in  the  frozen  and  thawed  condi¬ 
tion  (see  below).  Then,  TRANSFORM  reviews 
the  initial  part  of  the  FROST  output  file  to  deter¬ 
mine  the  material  type  identifier  for  each  element. 

From  the  incremental  (daily)  listings  of  the 
FROST  output  file,  TRANSFORM  next  reads  the 
following  conditions  for  each  element;  tempera¬ 
ture,  water  content,  ice  content,  material  density, 
and  porosity.  A  modulus  value  is  then  calculated 
for  each  element. 

Moduli  of  surface  paving  materials  were  calcu¬ 
lated  as  follows.  PCC  concrete  had  the  modulus 


set  at  a  constant  value  of  3.4  x  10^  kPa  (5,000,000 
Ib/in.^).  In  the  initial  two  simulation  series,  the 
resilient  modulus,  ,  of  asphalt  concrete  layers 
were  calculated  by  the  equation  (Schmidt  1975): 

M,(lb/in.2)  =  10(6-2«5-1.931xl0-2r 

-3.280x10*^  -1 ,888x10“’  T’ 

+1.175xl0-''r‘'+1.502xl0-*T’ 
-2.022xl0-^°r®) 

where  T  is  the  temperature  of  the  element  (°C). 
For  pavement  temperatures  greater  than  50°C,  the 
asphalt  modulus  was  set  to  1.7  x  10^  kPa  (25,000 
Ib/in.^);  for  temperatures  less  than  -29°C,  it  was 
set  to  3.3  X  10^  kPa  (4,840,000  Ib/in.^).  In  a  third 
series  of  simulations,  a  second  model  was  used 
for  predicting  asphalt  moduli  when  the  tempera¬ 
ture  was  above  1°C  (Ullidtz  1987): 

M,  (lb/in.2)  =  [15,000 

-7900  log  (T)]  X  145.04.  (9) 

At  below  1°C,  the  Schmidt  (1975)  relationship 
was  used.  Predictions  of  the  two  asphalt  models 
are  compared  in  Figure  4. 

In  the  layers  that  represent  an  unstabilized  base 
course,  subbase  or  subgrade  material,  TRANS¬ 
FORM  calculates  the  modulus  using  regression 
equations  developed  from  results  of  laboratory 
resilient  modulus  testing  conducted  on  frozen  and 
thawed  soil  samples  (Berg  et  al.  1996).  Each  soil 
element  is  first  classified  using  the  criteria  in  Table 
1  to  determine  which  type  of  modulus  equation  is 
appropriate.  The  modulus  is  then  calculated  using 
one  of  the  equations  shown  in  general  form  in 
Table  2.  The  equations  relate  the  frozen  resilient 
modulus  to  temperature  (through  unfrozen  water 
content),  and  the  unfrozen  resilient  modulus  to 
degree  of  saturation,  stress  condition,  and  to  den¬ 
sity  for  the  class  6  special  base  course  and  the 
1206  subgrade.  The  specific  parameters  for  the 
Mn/ROAD  materials  will  be  discussed  in  a  later 
section. 

The  following  parameters,  which  are  specific 
to  each  material,  are  used  in  the  regression  equa¬ 
tions  (Table  2):  1)  FCl  and  FC2,  regression  coef¬ 
ficients  for  the  resilient  modulus  in  the  frozen 
condition,  2)  TCI,  TC2,  TC3,  and  TC4,  regres- 


6 


(kPa)  (Ib/in^) 


Figure  4.  Comparison  of  predicted  values  from  the  Schmidt  and  Ullidtz 
asphalt  modulus  equations. 


Table  1.  Classification  of  soil  for  corresponding  modulus  equation. 

Classification _ Definition _ 

1  Frozen:  the  ice  content  is  greater  than  or  equal  to  0.005,  the  temperature  not  equal  to  0°C,  and  the 

pore  pressure  less  than  0. 

2  Partially  thawed:  the  ice  content  is  greater  than  or  equal  to  0.005  and  the  pore  pressure  is  greater 

than  or  equal  to  0. 

3  Thawed,  with  negative  pore  pressure:  the  ice  content  is  less  than  0.005  and  the  pore  pressure  is 

less  than  0. 

4  Thawed,  with  positive  pore  pressure:  the  ice  content  is  less  than  0.005  and  the  pore  pressure  is 

greater  than  0. 

5  Recently  frozen  layer:  the  temperature  is  less  than  0°C  and  was  recently  thawed.  This  layer 

retains  its  modulus  prior  to  freezing  until  the  frozen  equation  predicts  a  higher  modulus  or  until 
the  temperature  rises. 

6  Recovering  thawed  layer:  for  120  days  subsequent  to  having  been  frozen,  a  thawed  layer  has  its 

modulus  reduced  according  to  a  ratio  that  diminishes  over  time. 


sion  coefficients  for  the  resilient  modulus  in  the 
thawed  condition,  3)  /(CT)^^,  an  approximation  of 
the  mean  stress  condition,  and  4)  a  and  p,  unfro¬ 
zen  water  content  constants. 

Several  forms  of  the  unfrozen  water  content 
function  were  investigated  for  predicting  the  fro¬ 
zen  moduli  of  the  Mn/ROAD  materials,  as  re¬ 
ported  in  Berg  et  al.  (1996).  The  form  used  for  the 
Phase  1  simulation  series  of  this  study  is  the  nor¬ 
malized  gravimetric  unfrozen  water  content,  cal¬ 
culated  as  shown  in  Table  2.  A  second,  normal¬ 
ized  volumetric  form  of  the  unfrozen  water  content 
function  was  used  in  Phases  2  and  3  simulation 


series.  It  is  calculated  by  multiplying  the  gravi¬ 
metric  unfrozen  water  content  by  the  dry  density 
in  Mg/m^. 

TRANSFORM  uses  an  approximation  of  the 
/(a)  stress  term  in  its  equations  to  predict  moduli 
of  unfrozen  soils.  Since  this  term  is  constant  for 
each  soil,  the  unfrozen  moduli  calculated  by 
TRANSFORM  are  not  dependent  on  stress.  The 
approximation  used  was  the  mean  stress  condi¬ 
tion  applied  during  the  laboratory  resilient  modu¬ 
lus  testing. 

Because  some  materials  have  a  significant  dis¬ 
continuity  between  the  minimum  frozen  and 
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Table  2.  General  form  of  resilient  modulus  (Af^)  equations  used  in 
TRANSFORM. 


Condition 

Equation 

Frozen 

M,(lb/in.2)=FCl/(Wu_g)^<^2 

Partially  thawed 

M,(lb/in.2)  =  TC\  (100)^^V(Y)^‘^'‘/(<7)m^^ 

Thawed,  negative  pore  pressure 

M,(lb/in.2)  = 

Thawed,  positive  pore  pressure 

M;.(lb/in.2)  =  7’Cl(100)^‘^V(Y)^‘^‘*/(<^)m‘^^ 

Notes: 

FCl,  FC2  =  regression  coefficients  for  the  frozen  condition 


/K-g)  = 


a,  3  = 
7  = 
T  = 

TC1-TC4  = 


grav.  unfrozen  water  content  (decimal)  = 
unit  water  content  (1.0) 
unfrozen  water  content  constants 
temperature,  °C 
1.0°C 

regression  coefficients  for  the  thawed  condition 


T_ 

To 


f(S)  =  S/S^ 

S  =  degree  of  saturation  (%) 

S^=  1.0% 

/(Y)=  Ya/Yo 

Yj  =  dry  density  (Mg/m^) 

Y^  =  1.0  Mg/m^ 
a  =  stress  (Ib/in.^) 

/  (a)j^  =  a  mean  value  of  either  /j(g)  = 

OTf2(o)  =  ('^2^'^oc^^^o 
or/3(<^)  = 

a  =  1.0  lb/in-2 

o 

7-  =  bulk  stress  (lb/in.2)  =  30^  +  0^ 

—  2nd  stress  invariant  (Ib/in.^)  =  3a^  +  2 
=  octahedral  shear  stress  (Ib/in.^)  (-J2  /  3) 


thawed  predicted  moduli,  algorithms  were  added 
to  TRANSFORM  for  providing  a  smoother  tran¬ 
sition  in  calculated  moduli  between  the  frozen 
and  thawed  condition.  When  the  soil  goes  from  an 
unfrozen  to  a  frozen  state,  the  modulus  is  held 
constant  at  its  unfrozen  value  until  either  of  two 
changes  occur:  1)  the  soil  becomes  cold  enough 
that  the  modulus  as  calculated  by  the  frozen  equa¬ 
tion  exceeds  the  pre-frozen  modulus,  or  2)  the  soil 
warms.  When  the  soil  warms,  the  modulus  is  set 
at  the  value  calculated  by  the  frozen  equation  and 
continues  to  follow  the  frozen  curve  until  the  soil 


is  completely  thawed.  After  complete  thawing, 
the  soil  undergoes  a  120-day  “recovery  period.” 
During  this  time,  the  modulus  value  increases  ac¬ 
cording  to  a  power  relationship  with  about  one- 
half  of  the  recovery  taking  place  in  the  first  25 
days  about  75%  recovery  in  75  days  and  100% 
recovery  in  120  days. 

After  calculating  the  predicted  modulus  of  an 
element,  TRANSFORM  also  assigns  it  a  Poisson’s 
ratio  using  the  criteria  shown  in  Table  3. 

After  a  modulus  value  has  been  assigned  to  all 
elements  from  a  particular  time  step,  adjacent  ele- 
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Table  3,  Poisson’s  ratio  for  the  material  layers. 


Poisson ’s  ratio 

Condition 

Asphalt 

0.30 

Layer  temperature  is  less  than  -2.0°C 

0.35 

Layer  temperature  is  greater  than  or  equal  to  -2.0°C  and  less  than  or 
equal  to  1°C 

0.40 

Layer  temperature  is  greater  than  1°C  and  less  than  or  equal  to  8°C 

0.45 

Layer  temperature  is  greater  than  8°C 

Concrete 

0.15 

Constant  for  all  conditions 

Soil 

0.33 

Thawed,  volumetric  ice  content  is  less  than  0.005 

0.35 

Frozen,  volumetric  ice  content  is  greater  than  or  equal  to  0.005 

ments  within  the  same  material  type  are  combined 
into  a  single  layer  if  the  modulus  of  the  deeper 
element  is  less  than  ±20%  different  from  the  modu¬ 
lus  of  the  upper  sublayer.  A  “weighted  average” 
modulus  of  the  two  elements  is  then  determined, 
with  the  weighting  based  on  their  relative  lengths. 
The  modulus  of  the  next  lower  finite  element  is 
then  compared  with  the  modulus  of  the  upper  ele¬ 
ment.  The  checking  and  combining  process  con¬ 
tinues  until  an  element  modulus  is  outside  of  the 
20%  limitation  or  if  a  layer  of  a  different  material 
or  a  different  frozen/unfrozen  state  is  encountered. 
In  this  manner,  a  particular  material  layer  in  the 
pavement  profile  may  be  divided  into  several 
sublayers.  During  the  process  of  combining  ele¬ 
ments  with  similar  modulus  values,  the  thickness 
of  each  sublayer  is  also  determined,  as  well  as  a 
weighted  average  of  its  other  properties  such  as 
temperature,  density  and  Poisson’s  ratio. 

Typically,  the  pavement  profile  that  was  di¬ 
vided  into  99  finite  elements  for  the  FROST  pro¬ 
gram  is  combined  by  TRANSFORM  into  5  to  20 
sublayers  with  similar  resilient  modulus  values. 
During  the  winter  and  spring  a  larger  number  of 
sublayers  is  more  prevalent  than  in  the  summer 
months.  TRANSFORM  creates  an  additional  “in¬ 
finite”  layer  beneath  the  modeled  column  for  pass¬ 
ing  to  NELAPAV,  which  has  its  properties  set  the 
same  as  those  for  the  bottom  modeled  sublayer. 

Additional  items  must  be  input  to  (or  generated 
from)  TRANSFORM  in  order  for  it  to  create  in¬ 
put  files  for  NELAPAV,  They  include  1)  loading 
information  such  as  the  total  load  on  a  specified 
loaded  radius,  or  load  pressure,  2)  location  and 


depth  information  related  to  the  points  where  the 
stresses  and  strains  are  to  be  computed,  and  3)  a 
model  number  telling  NELAPAV  which  form  of 
the  modulus  equation  to  use  for  each  material. 

In  the  program,  a  4082-kg  (9000-lb)  load  was 
applied  to  a  radius  of  15.0  cm  (5.91  in.),  which 
approximates  the  area  of  a  standard  set  of  dual 
wheels  or  a  falling  weight  deflectometer  (FWD) 
testing  plate.  In  all  cases,  stresses  and  strains  were 
computed  beneath  the  center  of  the  load.  For  flex¬ 
ible  pavements,  the  stress  state  at  two  points  were 
analyzed:  at  the  bottom  of  the  pavement  layer, 
and  at  the  top  of  the  subgrade.  For  rigid  pave¬ 
ments,  stress  was  computed  only  at  the  bottom  of 
the  pavement.  In  all  cases  the  point  of  computa¬ 
tion  was  0.01  in.  from  the  interface  between  mate¬ 
rials. 

NELAPAV 

NELAPAV  is  an  acronym  for  Nonlinear  Elas¬ 
tic  Layer  Analysis  for  PAVements.  It  computes 
stresses,  strains,  and  displacements  at  any  point  in 
an  n-layered  pavement  system.  The  mainframe 
computer  version  of  the  program  was  developed 
by  Lynne  Irwin  of  Cornell  University  and  Gregor 
Fellers  of  CRREL  in  1980.  The  microcomputer 
version  was  developed  by  Irwin  and  Daniel  Speck 
at  Cornell  University  in  1984  and  1985  (Irwin  and 
Speck  1986).  The  program  is  an  adaptation  of  the 
Chevron  Layered  Elastic  Systems  program 
(CHEVLAY). 

Irwin  and  Speck  (1986)  describe  the  computa¬ 
tional  approach  used  by  NELAPAV  and  define 
the  following  terms.  The  term  state  of  a  point  in  a 
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layered  system  is  used  to  refer  to  the  set  of 
stresses,  strains,  and  displacements  occurring 
at  that  point  for  a  particular  loading  condition. 

A  set  of  layer  moduli  is  said  to  be 

compatible  at  radius  r  with  a  given  set  of  non¬ 
linear  models  (f if  each  M,  agrees 
(within  a  specified  tolerance)  with  the  modu¬ 
lus  determined  by  the  model  /;  evaluated  at 
the  middepth  of  layer  i  and  at  radius  r. 

The  fundamental  assumption  of  NELAPAV 
may  be  stated  as  follows:  “The  state  of  a  point 
p  =  (r,z)  in  a  layered  system  is  primarily  deter¬ 
mined  by  the  modulus  of  elasticity  of  each  of 
the  layers  in  the  system  at  the  radius,  r,  of  the 
point.”  Once  this  assumption  has  been  made, 
the  state  of  any  point  in  the  layered  system 
may  be  determined  by  computing  a  set  of  com¬ 
patible  moduli  at  the  radius  of  the  point  and 
then  evaluating  the  resulting  layered  system  at  the 
point,  holding  the  moduli  constant. 

NELAPAV  operates  as  a  “front-end”  to 
CHEVLAY  in  that  it  calls  the  basic  CHEVLAY 


(kPa)  (lb/in.2) 


Figure  5.  Stress  dependence  of  resilient  modulus  in  1206 
clay  subgrade  material  (from  Berg  et  al.  1996). 


routines  iteratively  in  order  to  come  up  with  a  set 
of  compatible  moduli  for  the  layered  system  at  the 
radius  of  each  of  the  desired  calculation  points. 
The  tolerances  for  compatibility  used  in  this  study 


Table  4.  Models  currently  available  in  NELAPAV  (from  Yang  1988). 

Model 


Name  _ Specification 


0 

Linear 

=  constant 

1 

Bulk  stress 

Mr  =  kiO'^^ 

2 

Deviator  stress 

(k2  +  k^(ki-C5^)  c^<ki 

M  =  \ 

3 

Second  stress  invariant 

4 

Octahedral  shear  stress 

Mf  =  ^ 

5 

Vertical  stress 

6 

Major  principal  stress 

Mr  =  A:iai*2 

7 

First  stress  invariant 
octahedral  shear  stress 
and  anisotropic 
consolidation  ratio 

Notes: 

0  =  bulk  stress 

“  constants 
=  deviator  stress 
a  =  vertical  stress 

V 

Oj  =  major  principal  stress 

=  first  stress  invariant  due  to  overburden  only 
J,  =  first  stress  invariant  due  to  overburden  and  load 

ip  ... 

=  anisotropic  consolidation  ratio 
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were  a  2.0%  maximum  change  in  from 
one  iteration  to  the  next  and  a  total  of  10 
iterations. 

NELAPAV  allows  the  use  of  nonlinear 
(i.e.,  stress-dependent)  modulus  values  in 
the  analysis.  Modulus  values  for  thawing 
and  unfrozen  fine-grained  soils  are  highly 
nonlinear  as  illustrated  in  Figure  5.  Table 
4,  from  Yang  (1988),  illustrates  the  vari¬ 
ous  types  of  linear  and  nonlinear  models 
currently  available  for  use  in  NELAPAV ; 
however,  models  2  and  7  are  not  currently 
incorporated  in  the  rest  of  the  CRREL 
design  procedure.  In  the  CRREL  version 
of  NELAPAV,  model  1  has  been  changed 
to  the  semi-log  form: 


Table  5.  Format  of  a  NELAPAV  input  file  (from  Irwin  and 
Speck  1986). 


Subject 

Header 


Layered  system 


Tolerance 


Calculation  points 


input  parameters 
problem  description 
input  units 

load,  pressure,  load  radius 
n 

M ^  E'j  seed  Vj  ^2,1  ^3,1  ^4,i 

£*2  seed  V2  Y2  ^0  2  ^2  2  ^2  2  ^3  2  ^4  2  h 


£  seed  v„  y„  n  n  ^ 

n  n  n  *n  0,n  n  l,n  2,n  3,  n  4,n  n 

maximum  iterations,  tolerance 


where  and  ^2  are  constants  and  0  is  bulk 
stress.  For  this  study,  we  utilized  models 
0,  1,3,  and  4,  the  specifics  of  which  will 
be  described  in  a  later  section.  - 

One  shortcoming  of  NELAPAV  is  that  Note 

it  allows  only  one  circular  load  to  be  ap¬ 
plied  at  the  surface.  This  is  not  a  major 
problem  in  this  analysis  because  we  are 
dealing  with  roadway  pavements  rather 
than  airport  pavements,  which  experience 
much  more  complex  tire  configurations. 

Table  5  is  a  brief  listing  of  the  informa¬ 
tion  included  in  a  NELAPAV  input  file  ^ 

(Irwin  and  Speck  1986).  It  may  contain 
up  to  25  layers,  with  the  following  values  r 

for  each:  model  number  (from  Table  4), 
seed  modulus,  Poisson’s  ratio,  total  den¬ 
sity,  lateral  earth  pressure  coefficient,  thickness, 
and  the  constants  required  for  the  equations  from 
Table  4. 

NELAPAV  requires  a  seed  modulus  to  begin 
its  predictions  of  modulus  for  the  various  material 
layers.  The  value  used  for  the  seed  modulus  is  the 
resilient  modulus  calculated  in  and  passed  from 
TRANSFORM.  The  constants  required  for  com¬ 
putations  are  also  passed  from  TRANSFORM. 
The  k2  value  is  a  constant  assigned  on  the  basis  of 
material  type.  The  constant  varies  with  the  mois¬ 
ture/density  level  of  the  layer,  and  is  calculated  by 
multiplying  together  all  the  nonstress  terms  in  the 
predictive  equation  for  the  material  (Table  2). 


Notes: 

n  =  number  of  layers 
Supplied  for  each  layer 

M  -  model  number  (see  Table  4) 

£  seed  =  seed  modulus  (Ib/in.^) 
v  =  Poisson’s  ratio 
y  =  density  (Ib/ft^) 

Kq  =  lateral  earth  pressure  coefficient 
h  =  thickness  (in.) 

ky.  k^  =  constants  to  be  used  in  models  (see  Table  4) 
t  =  mean  temperature  (°C) 

r  ,  z  =  radius,  depth  at  which  to  make  calculations  (in.) 


During  this  study,  ATq,  the  coefficient  of  lateral 
earth  pressure,  was  assigned  to  the  layers  as  fol¬ 
lows:  a  value  of  1.5  was  used  for  paving  materials 
and  frozen  soil;  a  value  of  1.0  was  used  for  all 
unfrozen  soil  layers.  We  now  feel  that  this  system 
is  too  simplified,  and  are  modifying  the  section  of 
TRANSFORM  that  assigns  Kq. 

For  each  layer,  NELAPAV  expects  to  receive  a 
value  for  the  total  density,  which  is  adjusted  to  a 
buoyed  density  below  the  water  table.  In  this  study, 
however,  we  used  dry  density  instead  of  total  den¬ 
sity,  and  neglected  to  adjust  to  a  buoyed  density 
below  the  water  table.  A  short  sensitivity  study, 
conducted  when  these  problems  were  discovered. 
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Table  6.  Cumulative  damage  models  used. 

A.  Flexible  Pavement  Horizontal  Strain  Criteria 

1)  The  Asphalt  Institute  (MS-1,  1982): 

iVa  =18.4  C  (4.325x10-^) 

where  =  number  of  load  applications  to  45%  cracking 

C  =  a  function  of  the  volume  of  the  voids  and  the  volume  of  asphalt,  10^ 
z  =  4.84  [(V^A^^  +  V^)- 0.69] 

=  volume  of  the  asphalt,  percent  (11%) 

V  =  volume  of  the  voids,  percent  (5%) 

Ej  =  tensile  strain  at  the  bottom  of  the  asphalt  layer,  in./in. 

=  modulus  of  the  asphalt  layer,  Ib/in.^ 

2)  Witczak  (1972)  (i.e..  Asphalt  Institute  MS-1 1); 

Aa  =  aM‘'(l/e,)" 

where  a  =  1.86351  x  10”'^ 
b  =  1.01996 
c  =  4.995 
d  =  1.45 

q  =  pavement  temperature,  °F 

3)  The  Corps  of  Engineers  (U.S.  Army  1988): 

^  jQ(2.68-51og|eJ-2.661og£:j 

4)  Coetzee  and  Connor  (1990): 

where  when  E^  ^  1,500,000  Ib/in.^:  a,  b,  c  =  3.364  x  10^,  —7.370,  —4.470 

and  when  E^  <  1,500,000  Ib/in.^:  a,  b,  c  =  6.565  x  10®,  -5.764,  -3.640 

B.  Flexible  Pavement  Subgrade  Strain  Criteria 

1)  The  Asphalt  Institute  (1982): 

_2Q[l/waog/-Iogev)] 

where  =  allowable  traffic  based  on  subgrade  strain 
m  =  a  constant  (0.25) 

I  =  a  constant  (2.8  x  10“^) 

=  vertical  strain  at  the  top  of  the  subgrade,  in./in. 

2)  The  Corps  of  Engineers  (U.S.  Army  1987): 
iVs  =10000  (A/E^f 

where  A  =  0.000247  +  0.00245  log 
B  =  0.0658 

=  subgrade  resilient  modulus,  IhUn? 


12 


Table  6  (cont’d). 

C.  Rigid  Pavement  Horizontal  Stress  Criteria 

1)  The  Corps  of  Engineers  (U.S.  1990): 


_  2Q  [(df-adon)/bdon] 


where:  =  allowable  traffic  based  on  horizontal  stress 

df  =  RJo, 

''  con  h 

R  -  flexural  strength  of  the  concrete,  Ib/in/ 

con  ®  ^ 

Ojj  =  horizontal  stress  at  the  base  of  the  concrete,  Ib/in/ 
adon  =  0.2967  +  0.002267  SCI 
bdon  =  0.3881  +0.000039  SCI 


SCI  =  surface  condition  index  of  the  pavement  when  failed 


showed  that  this  did  not  significantly  affect  the 
resultant  strain  calculations. 

The  solution  computed  by  NELAPAV  is  an 
approximation  of  the  exact  solution.  In  reality,  the 
stress  state  changes  from  point  to  point.  There¬ 
fore,  the  modulus  of  a  nonlinear  material  varies 
both  vertically  and  horizontally.  While  NELAPAV 
recomputes  the  set  of  compatible  moduli  to  deter¬ 
mine  the  states  of  points  at  different  radii,  it  is 
bound  by  the  assumption  that  the  moduli  are  con¬ 
stant  everywhere  in  the  layers.  An  exact  theory 
for  nonlinear  materials  would  allow  the  modulus 
to  vary  horizontally  within  the  layer  in  accordance 
with  the  nonlinear  model. 

Output  from  NELAPAV  is  one  file  for  each 
day.  The  files  include:  1)  a  repeat  of  the  input 
information,  2)  compatible  moduli  of  the  layers 
resulting  from  calculations,  and  3)  stress  condi¬ 
tions  for  all  points  specified.  For  this  study,  the 
points  specified  were  located  0.01  in.  above  the 
bottom  of  the  asphalt  or  PCC  and  0.01  in.  below 
the  top  of  the  subgrade. 

CUMDAM 

The  program  CUMDAM  calculates  cumulative 
damage  to  the  pavement  structure,  and  was  devel¬ 
oped  at  CRREL.  No  report  has  been  prepared  that 
discusses  its  function  and  operation.  In  general 
form,  the  procedure  used  for  CUMDAM ’s  calcu¬ 
lations  is  the  linear  summation  of  cycle  ratios, 
referred  to  as  Miner’s  mle,  which  may  be  stated  as: 


where  nj  =  number  of  applications  at  strain  level  i 
N■^  =  number  of  applications  to  cause  fail¬ 
ure  at  strain  level  i,  based  on  damage 
model  predictions 
D  =  total  cumulative  damage. 

In  this  relation,  failure  can  occur  when  D  equals 
or  exceeds  1 .0.  Thus,  for  a  section  to  last  its  de¬ 
sign  life,  the  value  of  D  should  not  accumulate  to 
1.0  until  the  design  period  expires. 

The  value  relates  to  the  design  traffic,  or 
applications,  in  8165-kg  (18,000- lb)  equivalent 
standard  axle  loadings  (ESALs).  For  the  first  two 
simulation  series  in  this  study,  anticipated  traffic 
was  considered  to  be  2,815,000  ESALs  during  a 
5-year  period.  Since  our  incrementation  was  on  a 
daily  basis,  we  applied  this  as  a  constant  loading 
of  1542  ESALs  per  day.  In  the  final  simulation 
series  the  anticipated  traffic  was  revised  upwards 
to  3,300,00  ESALs  over  five  years,  applied  at  a 
rate  of  1808  ESALs  per  day. 

CUMDAM  includes  several  damage  models 
previously  developed  by  others  for  determining 
N■^  the  number  of  applications  to  failure  at  a  par¬ 
ticular  strain/stress  condition  (Table  6).  Of  the 
damage  models  for  flexible  pavements,  four  are 
based  on  horizontal  strain  at  the  bottom  of  the 
asphalt  layer,  and  relate  to  damage  effects  that 
result  in  pavement  cracking.  These  models  were 
developed  by  the  Asphalt  Institute  ( 1 982),  Witczak 
(1972),  the  Corps  of  Engineers  (U.S.  Army  1988), 
and  Coetzee  and  Connor  ( 1 990).  Three  other  dam¬ 
age  models  for  flexible  pavements  are  based  on 


13 


the  vertical  strain  at  the  top  of  the  subgrade,  and 
these  relate  to  rutting  damage  in  the  pavement. 
They  were  developed  by  the  Asphalt  Institute 
(1982),  the  Federal  Aviation  Administration  (Bush 
1980),  and  the  Corps  of  Engineers  (U.S.  Army 
1987).  For  rigid  pavements,  CUMDAM  uses  the 
damage  model  developed  by  the  Corps  of  Engi¬ 
neers  (U.S.  Army  1990),  which  is  based  on  the 
horizontal  stress  at  the  base  of  the  PCC. 

The  program  assumes  that  all  applications  will 
affect  the  point  being  modeled.  That  is,  the  dam¬ 
age  is  not  reduced  according  to  a  pass-to-cover- 
age  algorithm  to  simulate  the  lateral  wander  of  the 
axle  within  the  travel  lane. 

It  is  recognized  that  some  of  the  equations  are 
being  applied  outside  of  the  original  assumptions 
used  in  their  development;  however,  they  are  rep¬ 
resentative  of  cumulative  damage  models  currently 
available  and  are  used  for  the  initial  analysis  until 
more  appropriate  equations  can  be  developed. 

The  CUMDAM  program  first  reads  the  strain/ 
stress  conditions  from  the  daily  NELAPAV  out¬ 


put  files  and  calculates  the  allowable  applications 
predicted  by  the  various  models  under  those  con¬ 
ditions.  It  then  divides  the  design  applications  by 
the  allowable  applications,  producing  a  daily  in¬ 
cremental  damage  value.  The  incremental  values 
are  summed  to  produce  a  cumulative  damage  value 
and  these  are  printed  to  individual  files  for  each 
model  type.  A  file  summarizing  the  daily  strain 
and  deflection  information  is  also  produced. 


MN/ROAD  PAVEMENT 
PERFORMANCE  STUDIES 

Performance  predictions  using  the  Mechanistic 
Pavement  Design  Procedure  were  conducted  in 
three  major  efforts: 

Phase  1,  conducted  in  the  spring  of  1991,  in¬ 
cluded  an  initial  simulation  series  that  modeled 
temperatures  from  a  year  close  to  the  mean  freez¬ 
ing  index.  These  boundary  conditions  were  ap¬ 
plied  to  eight  flexible  and  three  rigid  sections. 


Table  7.  Layer  composition  and  thicknesses  of  pavement  structure  in  test  sec 
tions  simulated. 


Test 

Section 

Layer  1 

Layer  2 

Layer  3 

Comp 

*  Thick  ^ 

Comp 

Thick 

Comp 

Thick 

Flexible 

5-yr  design 

F-1 

AC 

14.6  (5.75) 

CL4 

83.8  (33.0) 

F-2 

AC 

14.6  (5.75) 

CL6 

10.2  (4.0) 

CL4 

71.1  (28.0) 

F-3 

AC 

14.6  (5.75) 

CL5 

10.2  (4.0) 

CL3 

83.8  (33.0) 

F-4 

AC 

22.2  (8.75) 

10-yr  design 

F-14 

AC 

27.3  (10.75) 

F49 

AC 

19.7  (7.75) 

CL3 

71.1  (28.0) 

F-21 

AC 

19.7  (7.75) 

CL5 

58.4  (23.0) 

F-22 

AC 

19.7  (7.75) 

CL6 

45.7  (18.0) 

Rigid 

5-yr  design 

R-5 

PCC 

19.7  (7.5) 

CL4 

7.6  (3.0) 

CL3 

68.5  (27.0) 

R-6 

PCC 

19.7  (7.5) 

CL4 

12.7  (5.0) 

lO-yr  design 

R-ll 

PCC 

24.1  (9.5) 

CL5 

12.7  (5.0) 

*Composition: 

AC — asphalt  concrete;  PCC — portland  cement  concrete;  CL3 — class  3  special; 
CL4 — class  4  special;  CL5 — class  5  special;  CL6 — class  6  special 
^Thickness,  cm  (in.) 
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I  Asphalt  Concrete  Portland  Cement  Concrete  □  Class  3 


(cm)  (in.)  Class  4  Class  5  Class  6 


F-1  F-2  F-3  F-4  F-14  F-19  F-21  F-22 

5  Year  (ML  5-)  |  10  Year  (ML  10-) 

R-5  R-6  j  R-11 

5  Year  |  10  Year 

Flexible 

Rigid 

Figure  6.  Pavement  structure  ofMn/ROAD  test  sections  simulated. 


Phase  2,  an  effort  in  the  summer  of  1992,  had 
two  primary  objectives.  Phase  2A  included  three 
series  modeling  the  eight  flexible  sections  with 
the  mean  freeze  season  and  changing  the  method 
employed  to  calculate  the  asphalt  and  subgrade 
modulus.  Phase  2B  investigated  the  variability  in 
predictions  when  temperatures  from  freeze  sea¬ 
sons  with  maximum  and  minimum  freezing  indi¬ 
ces  are  applied  to  a  single  flexible  section. 

Phase  3,  an  effort  in  the  summer  of  1993,  ex¬ 
panded  the  investigation  of  the  effects  of  freeze 
season  characteristics.  This  series  modeled  21  dif¬ 
ferent  freeze  seasons  applied  to  one  full-depth  and 
one  conventional  flexible  section. 

Phase  1 

Pavement  sections 

The  Mn/ROAD  mainline  facility  has  23  test 
sections,  including  both  flexible  and  rigid  pave¬ 
ment  systems,  that  are  designed  to  fail  after  either 
5  or  10  years.  Included  in  the  test  matrix  are  pave¬ 
ment  systems  with  various  thicknesses  and  quali¬ 
ties  of  base  and  subbase  materials.  The  base/sub¬ 
base  materials  have  a  “special”  Mn/DOT  grading 
designation  for  the  test  facility  and  range  from 
class  3,  a  well-graded  sand,  through  classes  4  and 
5,  which  have  decreasing  percentages  of  fines,  to 
class  6,  a  well-graded  gravel  with  sand.  The  sub¬ 
grade  at  the  site  classifies  as  a  sandy  lean  clay,  a 


CL  in  the  Unified  Soil  Classification  System  emd 
an  A-6  in  the  AASHTO  Classification  System 
(Bigl  and  Berg  1996a).  Of  the  four  subgrade 
samples  tested  at  CRREL,  two  exhibited  high  heav¬ 
ing  in  the  frost  susceptibility  test,  while  the  other 
two  ranked  as  having  a  medium  heave  rate. 

The  mechanistic  procedure  was  applied  to  1 1 
test  sections,  which  were  distributed  among  the 
designs  as  follows:  four  5-year  flexible,  two  5- 
year  rigid,  four  10-year  flexible,  and  one  10-year 
rigid.  Figure  6  and  Table  7  show  the  pavement 
structures  of  the  test  sections  simulated.  Seven  of 
the  test  sections  had  variable  water  table  depths 
along  their  lengths,  and  these  were  simulated  at 
the  two  extreme  water  table  positions.  Table  8 
lists  the  test  sections  and  water  table  conditions 
simulated  as  well  as  the  nomenclature  used  to 
denote  the  various  cases.  The  case  nomenclature 
includes  the  Mn/DOT  test  section  number  and  the 
water  table  depth  (in  ft)  preceded  by  a  “w.”  For 
example,  flw9  refers  to  Mn/ROAD  test  section 
ML5-F-1  with  a  water  table  2.7  m  (9  ft)  below  the 
pavement  surface  and  rllw6  refers  to  section 
MLlO-R-11  with  a  water  table  1.8  m  (6  ft)  below 
the  pavement  surface.  All  but  one  case  involved 
the  high-heaving  subgrade  (sample  1206)  beneath 
the  pavement  structure.  The  final  simulation  in¬ 
cluded  the  lower-heaving  subgrade  (sample  1232) 
under  the  5 -year  full-depth  section  (ML5-F4), 
termed  case  f4w6ss,  or  “second  subgrade.” 
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Table  8.  Test  sections  and  conditions  analyzed. 


Water  table 


Mn/ROAD 
test  section 

Depth 

m(ft) 

Simulation 

case 

Other  variations 

FLEXIBLE 

5-year 

ML5-F-1 

2.7  (9) 

flw9 

ML5-F-2 

2.7  (9) 

f2w9 

ML5-F-2 

6.1  (20) 

f2w20 

ML5-F-3 

2.4  (8) 

f3w8 

ML5-F-3 

6.1  (20) 

f3w20 

ML5-F-4 

1.8  (6) 

f4w6 

ML5-F-4 

1.8  (6) 

f4w61d 

Low  density  subgrade 

ML5-F-4 

1.8  (6) 

f4w6ss 

No.  1232  subgrade 

10-year 

MLlO-F-14 

15.3  (50) 

fl4w50 

MLlO-F-19 

5.5  (18) 

fl9wl8 

MLlO-F-21 

5.5  (18) 

f21wl8 

MLlO-F-21 

13.7  (45) 

f21w45 

MLlO-F-22 

13.7  (45) 

f22w45 

RIGID 

5-year 

ML5-R-5 

1.8  (6) 

r5w6 

ML5-R-5 

3.7  (12) 

r5wl2 

ML5-R-6 

1.2  (4) 

r6w4 

ML5-R-6 

3.7  (12) 

r6wl2 

10-year 

MLlO-R-ll 

1.8  (6) 

rllw6 

MLlO-R-11 

3.7  (12) 

rllwl2 

Material  properties 

Material  properties  input  to  the  FROST  pro¬ 
gram  are  shown  in  Table  9.  The  soil  water  charac¬ 
teristics  are  based  on  a  best  fit  to  the  moisture 
retention  test  data  using  the  Gardner’s  eq  2.  Soil 
residual  moisture  contents  were  calculated  using 
eq  2  with  a  pressure  of  -800  cm  water.  The  soil 
permeability  characteristics  result  from  fitting  the 
Gardner’s  eq  3  to  unsaturated  hydraulic  conduc¬ 
tivity  test  results. 

As  part  of  this  study,  Mn/DOT  class  3  special, 
class  6  special,  and  two  subgrade  materials  had 
been  tested  to  determine  the  information  for  input 
to  the  model.  Results  of  the  soil  characterization 
tests  are  given  in  Bigl  and  Berg  (1996a).  The  class 
4  special  and  class  5  special  materials  were  not 
tested,  so  their  behavior  had  to  be  approximated 
using  data  from  previously  tested  materials  that 
most  closely  matched  their  specified  size  grada¬ 
tions.  A  subbase  from  taxiway  A  at  the  Albany, 


New  York,  airport  most  closely  matched  the  class 
4  special  subbase  specifications;  dense-graded 
stone,  which  had  been  tested  during  a  cooperative 
study  in  Winchendon,  Massachusetts,  most  closely 
matched  the  class  5  special  material.  Results  of 
the  characterization  tests  of  the  substitute  materi¬ 
als  are  also  given  in  Bigl  and  Berg  (1996a)  and  in 
Cole  etal.  (1986, 1987). 

The  densities  of  the  Mn/ROAD  materials,  which 
had  been  characterized  with  compaction  testing, 
were  set  at  the  optimum  values.  For  the  substitute 
materials,  the  densities  were  set  at  the  values  de¬ 
termined  from  in-place  samples  taken  from  the 
site.  To  demonstrate  the  influence  of  proper  com¬ 
paction,  one  simulation  (case  f4w61d)  involved 
setting  the  density  of  the  1206  subgrade  to  a  re¬ 
duced  value  less  than  optimum  (1.69  Mg/m^,  1 05.5 
Ib/ft^).  Optimum  density  was  set  at  1.89  Mg/m^ 
(117.9  Ib/ft^),  determined  with  a  compactive  ef¬ 
fort  of  2360  kJ/m3  (55,000  ft-lb/ft^). 
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Table  10.  Modulus  equations  used  for  Mn/ROAD  test  section  simulations. 


A.  Frozen  condition 

Std. 

Material 

Equation  (M^  in  Ib/in?)"^ 

n 

r2 

error 

Clay  subgrade  sample  1206  (565) 
Frozen 

Mr=l,049/(Wu_gr2-344 

207 

0.99 

0.275 

Frozen 

M,=l,052/(Wu_vr^-^^^ 

207 

0.99 

0.262 

Clay  subgrade  sample  1232  (566) 

Frozen 

Mr=846/(w„_g)-2*6l 

244 

0.98 

0.423 

Frozen 

M,=848/(w„_vr2'^” 

244 

0.98 

0.394 

Class  3  special  “stockpile” 

Frozen 

Mr=5,488/(Wu_grl<^’® 

210 

0.97 

0.507 

Frozen 

Mr=5,542/(w„_vrl-249 

186 

0.97 

0.467 

Class  4  special  (taxiway  A  subbase) 
Frozen 

Mr  =  l,813/K_gr‘''33 

85 

0.93 

0.885 

Frozen 

Mr=l,652/(wu_vr2-^'^ 

69 

0.91 

0.916 

Class  5  special  (dense  graded  stone) 

Frozen 

Mr=8,695/(Wu_g)-‘-28i4 

28 

0.95 

0.511 

Frozen 

M,=9,245/(w„_^r*'^39 

28 

0.97 

0.432 

Class  6  special  “stockpile” 

Frozen 

M,  =  19,427/(wu_g)-0-’95 

260 

0.98 

0.338 

Frozen 

Mr=19,505/(w„_^)-o«97 

260 

0.98 

0.341 

Notes: 

n  =  number  of  test  points 

=  coefficient  of  determination 
M  =  resilient  modulus 
AS)  =  s/s^ 

S  =  degree  of  saturation  (%) 

5^  =  1.0% 

/(y)  =  Yo 

=  dry  density  (Mg/m^) 

=  l.OMg/m^ 

f(w  )  =  w  /w^ 

u— u— g  o 

^  =  gravimetric  unfrozen  water  content 

w  =  unit  water  content  (1.0) 
o 

w  =  volumetric  unfrozen  water  content 


a  =  stress  (Ib/in.^) 

Aia)  = 
f^(o)  = 

=  l.Olb/in.2 

/j  =  bulk  stress  (Ib/in.^) 

7j  =  3  a^+  Gj 

=  2nd  stress  invariant  (Ib/in.^) 
h  ~  303^  +  203©^ 

=  octahedral  shear  stress  (Ib/in.^) 

C  - 


*Output  from  equations  can  be  converted  to  kilopascals  by  multiplying  by  6.895. 


Saturated  hydraulic  conductivities  were  set  at 
the  laboratory-determined  value,  except  in  the  case 
of  the  1206  clay  subgrade.  In  that  case,  we  multi¬ 
plied  the  laboratory  value  by  10  to  approximate 
the  hydraulic  behavior  after  formation  of  frost- 
induced  shrinkage  cracks. 


The  modulus  equations  used  in  TRANSFORM 
for  the  Phase  1  modeling  are  the  normalized  gravi¬ 
metric  form  shown  in  Table  10a  and  the  unfrozen 
form  shown  in  Table  10b.  The  associated  unfro¬ 
zen  water  constants  utilized  in  the  frozen  form  of 
the  equations  are  given  in  Table  11.  Development 
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Resilient  Modulus 


Table  10  (cont’d). 


B.  Unfrozen  condition 


Material 

Equation  (M^  in  Ib/in.^)* 

n 

r2 

Std. 

error 

Clay  subgrade  sample  1206  (565) 
Never  Frozen 

M,  =  1,597, 000/(5)-2-63/(y)14.42j^((j)-O.257 

655 

0.82 

0,251 

Clay  subgrade  sample  1232  (566) 
Never  Frozen 

M,  =  1.518  X  103°/(5)-’3®5/3((j)-0.272 

451 

0.95 

0.328 

Class  3  special  ^^stockpile” 

Thawed 

Mr  =283, 300/(5)-* 

408 

0.86 

0.520 

Class  4  special  (taxiway  A  subbase) 
Thawed 

M,  =  8.946  X  10V('5)-^°2®/2(o)°-292 

149 

0.86 

0.168 

Class  5  special  (dense  graded  stone) 
Thawed 

Mr  =  382,400/(S)-°*’59^2(a)0*®^ 

64 

0.77 

0.164 

Class  6  special  ‘‘stockpile” 

Thawed 

Mr  =  1,391/(5)-®-50’7/(Y)'*-°^/i(o)°®°* 

492 

0.79 

0.232 

Thawed 

Mr  =  5, 257/(5)-®-^*®  /(Y)'*®^  ^  0.0193/i  (o) 

492 

0.76 

0.249 

(kPa)  (Ib/in.^) 


Figure  7.  Predicted  moduli  ofMn/ROAD  materials. 


of  these  equations  is  discussed  in  Berg  et  al.  (1996). 
Figure  7  shows  the  predicted  moduli  for  the  Mn/ 
ROAD  materials  using  mean  values  for  stress  and 
density,  where  applicable. 

It  should  be  noted  that  the  thawed  modulus 
equation  used  for  class  6  special  material  is  an 
exception  to  general  form  given  in  Table  2.  In  this 
case,  we  used  a  semi-log  form  (eq  10),  rather  than 
the  log-log  form,  which  eliminated  problems  that 


occurred  from  negative  stresses  generated  by  the 
model. 

It  should  also  be  noted  that  after  the  Phase  1 
simulations  were  completed,  we  discovered  that 
resilient  modulus  tests  on  the  1206  subgrade  in 
the  unfrozen  condition  were  probably  in  error  due 
to  a  miscalibrated  testing  system.  As  a  result,  the 
unfrozen  subgrade  moduli  predicted  in  Phase  1 
modeling  are  likely  to  be  substantially  higher  than 
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Table  11.  Constants  for  equations  to  de¬ 
termine  gravimetric  unfrozen  moisture 
content,  H'u-g* 


Equation  constants'^ 
Soil _ a _ p 

Subgrade 


1206 

11.085 

-0.274 

1232 

8.121 

-0.303 

Class  3  special 

stockpile 

1.497 

-0.709 

Class  4  special 

Taxiway  A 

3.0 

-0.25^ 

Class  5  special 

Dense  stone 

2.0 

-0.40' 

Class  6  special 
Stockpile 

0.567 

-1.115 

*  Wu-g,  decimal  form  =  ^ 

^  Values  for  these  materials  are  estimated 


To 


exist  in  the  field,  resulting  in  less  damage  than 
would  have  been  obtained  with  more  reasonable 
moduli. 

The  model  numbers  for  the  program  NELAPAV 
(Table  4)  were  set  as  follows.  Model  0,  the  linear 
model,  was  used  at  all  times  for  paving  materials 
and  for  soils  in  the  frozen  state.  With  model  0,  the 
seed  modulus  is  accepted  as  the  modulus  of  the 
material  layer.  Model  numbers  were  assigned  to 
unfrozen  soils  as  follows:  Model  1 — class  6;  Model 
3 — class  3  and  class  4;  Model  4 — both  clay 
subgrades. 

In  all  of  the  cumulative  damage  studies  on  the 
flexible  pavements,  it  was  assumed  that  the  pave¬ 
ment  properties  were  constant  from  one  pavement 
test  section  to  another.  Some  of  the  test  sections 
we  simulated  will  have  pavement  properties  dif¬ 
ferent  from  others  because  of  the  experimental 
design  of  the  paved  surface.  Some  variations  in 
the  asphalt  pavement  properties  could  be  consid¬ 
ered  in  subsequent  modeling  efforts,  but  we  main¬ 
tained  the  same  asphalt  pavement  properties  in  all 
of  these  simulations. 


(cm)  (in.) 


Figure  S.  Example  output  from  FROST  for  Mn/ROAD  test  sec¬ 
tion  ML5-F4  with  a  L8-m  (6-ft)  water  table.  Simulation  starts 
on  1  October. 


Results— flexible  sections 

Output  from  FROST  for  the  5-year 
full-depth  simulation  (f4w6)  is  shown 
in  Figure  8.  The  top  graph  shows  the 
daily  mean  air  temperature  for  the  pe¬ 
riod  from  1  October  1959  to  14  No¬ 
vember  1960.  The  center  portion  il¬ 
lustrates  the  predicted  frost  heave  and 
the  bottom  graph  contains  the  pre¬ 
dicted  frost  and  thaw  penetration  as 
functions  of  time.  Frost  output  graphs 
for  all  the  flexible  cases  are  compiled 
in  Appendix  A.  Table  12  contains 
maximum  frost  heave  and  maximum 
frost  penetration  depths  for  each  of 
the  cases  simulated  (Table  8).  Simu¬ 
lations  with  shallower  water  table 
depths  had  greater  amounts  of  heave 
and  less  frost  penetration  compared 
to  simulations  with  a  deeper  water 
table  location.  Sections  that  included 
the  substitute  for  the  class  4  special 
subbase  (taxiway  A  subbase)  had 
higher  amounts  of  heave  than  those 
with  other  combinations  of  subbase 
materials. 
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Table  12.  Maximum  frost  heave  and  frost 
penetration  in  modeled  simulations 


Frost 

Frost 

Cases 

heave 

penetration 

cm  (in,) 

cm  (in.) 

5- YEAR 

Flexible 


flw9 

0.25  (0.10) 

137.9  (54.3) 

f2w9 

1.19  (0.47) 

121.9  (48.0) 

f2w20 

0.33  (0.13) 

142.0  (55.9) 

f3w8 

0.08  (0.03) 

134.1  (52.8) 

f3w20 

0.00 

142.0  (55.9) 

f4w6 

1.37  (0.54) 

118.1  (46.5) 

f4w6  Id 

1.32  (0.52) 

118.1  (46.5) 

f4w6  ss 

0.03  (0.01) 

137.9  (54.3) 

Rigid 

r5w6 

1.04  (0.41) 

121.9  (48.0) 

r5wl2 

0.00 

137.9  (54.3) 

r6w4 

4.72  (1.86) 

103.9  (40.9) 

r6wl2 

0.00 

134.1  (52.8) 

10- YEAR 

Flexible 

fl4w50 

0.00 

134.1  (52.8) 

fl9wl8 

0.00 

142.0  (55.9) 

f21wl8 

0.00 

134.1  (52.8) 

f21w45 

0.00 

134.1  (52.8) 

f22w45 

0.00 

134.1  (52.8) 

Rigid 

rl  lw6 

1.42  (0.56) 

118.1  (46.5) 

rllwl2 

0.00 

134.1  (52.8) 

Id  =  low  density  (1.69  Mg/m^ 

or  105.5  lb/  ft3) 

ss  =  second  subgrade  (1232) 

Figure  9  presents  the  moduli  being  calculated 
by  TRANSFORM  and  passed  as  seed  moduli  to 
NELAPAV  for  the  f4w6  case.  The  plotted  moduli 
are  the  minimum  values  for  each  day,  with  the 
subgrade  (top)  being  located  within  0.3  m  (1  ft)  of 
the  asphalt  and  subgrade  (bottom)  being  the  rest 
of  the  modeled  section  (to  4.0  m  or  13.1  ft).  Note 
that  the  predicted  asphalt  moduli  during  the  sum¬ 
mer  months  are  smaller  than  those  of  the  subgrade, 
which  results  in  high  horizontal  strains  at  the  base 
of  the  asphalt  layer  during  the  summer  months.  In 
later  simulations  (Phase  2  and  3),  we  used  other 
models  to  determine  the  asphalt  modulus. 

Figure  1 0  shows  the  deflection  and  strains  com¬ 
puted  by  NELAPAV  for  the  f4w6  case.  The  hori¬ 
zontal  strain  at  the  base  of  the  asphalt  is  low 
through  the  winter  with  a  few  small  peaks  during 
short-term  thaw  periods,  then  exhibits  high  peaks 
during  spring  thaw  and  summer.  It  appears  to  be 
closely  tied  to  the  asphalt  temperature/modulus 
during  the  spring  thaw.  The  vertical  strain  at  the 
top  of  the  subgrade  is  also  low  during  the  winter 
with  three  sharp  peaks  during  thaw  events.  The 
vertical  strain  also  rises  during  spring  and  sum¬ 
mer,  closely  tracking  the  pavement  surface  tem¬ 
perature.  The  deflection  plotted  is  from  the  base 
of  the  asphalt,  but  is  assumed  to  be  similar  to  the 
surface  deflection.  It  also  shows  peaks  during  mid¬ 
winter  thaws,  and  a  maximum  peak  during  the 
spring  prior  to  drainage  of  excess  moisture.  Gen¬ 
erally,  deflections  in  the  summer  are  greater  than 
winter  values. 


(kPa)  (lb/in.2) 


Figure  9.  Seed  moduli 
output  by  TRANSFORM 
for  Mn/ROAD  test  section 
ML5-F4  with  a  1.8-m  (6- 
ft)  water  table.  Sub  grade 
(top)  is  within  the  upper 
0.3  m  (1  ft)  of  sub  grade 
beneath  the  asphalt.  Sub¬ 
grade  (bottom)  is  beneath 
that  layer  to  the  bottom 
of  the  modeled  section. 
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Figure  10.  Deflection  and  strain  calculated  by  NELAPAV 
for  Mn/ROAD  test  section  ML5-F4  with  a  1.8-m  (6-ft) 
water  table.  Horizontal  strain  is  at  the  base  of  the  pave¬ 
ment.  Vertical  strain  is  at  the  top  of  the  subgrade. 


Table  13  lists  the  applications  to  failure  pre¬ 
dicted  from  the  1-year  simulation  period  ap¬ 
plied  to  the  flexible  sections.  The  values  listed 
in  Table  13  were  determined  by  taking  the  re¬ 
ciprocal  of  the  one  year  (365 -day)  cumulative 
damage  value  produced  by  the  CUMDAM  pro¬ 
gram  and  multiplying  by  the  number  of  appli¬ 
cations  received  by  the  sections  during  that 
time  period  (562,830  per  year).  This  applica¬ 
tion  rate  was  based  on  an  initial  estimate  that 
the  traffic  across  the  test  sections  would  be 
2,815,000  ESALs  in  5  years,  or  5,630,000 
ESALs  in  10  years.  Note  that  the  values  listed 
in  Table  15  are  thousands  of  ESAL  applica¬ 
tions.  Subsequent  to  conducting  the  Phase  1 
series,  a  revised  estimate  is  that  the  5-year  sec¬ 
tions  will  receive  3,300,000  ESAL  applications 
in  five  years,  (i.e.,  are  designed  to  fail  after 
3,300,000  ESAL  applications);  the  10-year  sec¬ 
tions  will  fail  at  6,600,000  ESAL  applications. 

Individual  values  listed  in  Table  13  vary 
anywhere  from  very  much  below  to  very  much 
above  the  failure  values  described.  By  compar¬ 
ing  values  horizontally  across  the  table  for  any 
case,  it  can  be  seen  that  different  models  esti- 


Table  13.  Predicted  applications  to  failure  (x  1000)  from  Phase  1  simulation  series  of  flexible 
pavement  test  sections. 


Horizontal  criteria _  _ Vertical  criteria 


Cases 

Asp.  Inst. 
MS-1 

Asp.  Inst. 
MS-11 

Corps  of 
Engineers 

Coetzee/ 

Connor 

Asp.  Inst. 
MS-1 

Corps  of 
Engineers 

FAA 

5-YEAR 

F1W9 

86 

498 

88 

4,380 

136,942 

5,653 

11,879 

F2W9 

522 

1,244 

550 

42,703 

32,477 

1,345 

1,526 

F2W20 

111 

1,534 

764 

66,450 

155,909 

8,226 

9,855 

F3W8 

957 

1,706 

958 

92,419 

116,287 

2,653 

9,611 

F3W20 

1,068 

1,801 

1,055 

105,202 

222,462 

7,889 

20,632 

F4W6 

28,585 

75,751 

90,926 

>10'^ 

3,715 

15 

152 

F4W61d 

1,905 

20,209 

7,036 

>10'^ 

165 

>10’ 

>10’ 

F4W6SS 

72 

1,180 

96 

>10’ 

>10’ 

>10’ 

>10’ 

10-YEAR 

F14W50 

103,272 

485,198 

605,194 

>  10’ 

8,602 

10 

251 

F19W18 

396 

6,148 

966 

75,751 

2,962,263 

56,680 

457,585 

F21W18 

3,394 

30,456 

12,751 

2,010,107 

562,830 

84,382 

42,478 

F21W45 

3,404 

30,358 

12,794 

2,010,107 

611,772 

92,571 

48,188 

F22W45 

7,614 

57,315 

31,147 

8,040,429 

73,190 

22,415 

3,388 

Notes: 

Id  =  low  density  (1.69  Mg/m^  or  105.5  Ib/ft^) 
ss  =  second  subgrade  (1232) 

Traffic  simulated  at  rate  of  562,830  ESAL  applications/yr 
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Elapsed  Time  (doys) 


Figure  11.  Cumulative  damage  for  case  f4w6  with  opti¬ 
mum  density  1206  subgrade. 


Elapsed  Time  (doys) 


Figure  12.  Cumulative  damage  for  case  f2w9. 


Figure  13.  Cumulative  damage  for  case  f4w6  with  1232 
subgrade. 


mate  substantially  different  pavement  lives  for  the 
same  strain  conditions. 

A  close  review  of  the  values  in  Table  13  illus¬ 
trates  a  variety  of  interesting  occurrences; 

1.  Cases  with  Mn/ROAD  test  section  ML5-F-4 


that  employed  the  lower  density  subgrade 
(f4w61d)  or  the  “1232”  subgrade  (f4w6ss)  pre¬ 
dicted  the  pavement  to  fail  much  more  quickly 
than  when  the  higher  (optimum)  density  1206 
subgrade  was  used.  This  is  due  to  the  fact  that 
the  higher  density  “1206”  subgrade  used  in 
most  simulations  maintained  a  resilient  modu¬ 
lus  during  the  summer  months  in  excess  of 
276,000  kPa  (40,000  Ib/in.^),  but  the  modulus 
values  for  the  same  period  were  about  96,000 
kPa  (14,000  Ib/in.^)  for  f4w61d  and  about  6900 
kPa  (1000  Ib/in.^)  for  f4w6ss.  The  “low  den¬ 
sity”  predicted  moduli  are  the  closest  to  the 
values  measured  during  the  summer  months 
with  an  FWD. 

2.  Variations  in  water  table  position  in  the 
same  test  section  result  in  different  predicted 
lifespans,  especially  if  the  water  table  is  less 
than  3  m  (10  ft)  deep  at  one  location  and  more 
than  6  m  (20  ft)  deep  at  another.  Both  the 
ML5-F-2  and  ML5-F-3  test  section  simula¬ 
tions  illustrate  a  shorter  predicted  life  with  a 
shallow  water  table  depth.  For  the  shallow 
water  table  conditions,  greater  availability  of 
water  results  in  higher  thaw  weakening  (re¬ 
duced  modulus  values)  in  the  spring. 

3.  Most  of  the  sections  are  predicted  to  fail 
due  to  cracking  of  the  pavement  caused  by 
tensile  stresses  at  the  bottom  of  the  pavement 
much  before  they  will  fail  due  to  rutting  caused 
by  excessive  deformation  of  the  subgrade.  It 
is  only  the  simulations  of  the  full-depth  sec¬ 
tions  (ML5-F-4,  case  f4w6,  and  MLlO-F-14, 
fl4w50)  that  indicate  rutting  failures  may  oc¬ 
cur  before  excessive  fatigue  cracking. 

Figures  11  to  13  illustrate  examples  of  the 
relation  between  accumulating  damage  and 
time.  Additional  plots  of  heave,  frost  penetra¬ 
tion,  and  cumulative  damage  for  all  the  flex¬ 
ible  sections  are  compiled  in  Appendix  A. 

Damage  related  to  the  horizontal  strain  cri¬ 
teria  has  different  patterns  with  time  for  dif¬ 
ferent  models.  Damage  predicted  by  the  MS-1 
(Asphalt  Institute)  model  has  some  increase 
during  winter  thaw  events  and  during  the  main 
spring  thaw,  but  a  majority  of  its  increase  is  in  the 
summer  period.  Damage  from  the  MS-11  (As¬ 
phalt  Institute)  model  has  a  consistent  pattern  of 
rising  during  winter  and  spring  thaw  periods,  and 
remaining  constant  during  the  summer.  The  pat- 
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Table  14.  Applications  to  failure  (x  1000)  from  simulations  of 
rigid  test  sections. 


RCON  =  650  RCON  =  500  RCON  =  650 


Case 

SCI  =  80 

SCI  =  80 

SCI  =  90 

5-YEAR 

R5W6 

52,601 

507 

45,244 

R5W12 

65,598 

605 

56,452 

R6W4 

157,655 

1,505 

135,622 

R6W12 

132,120 

1,701 

113,933 

10-YEAR 

R11W6 

>5.6  X  10* 

18,761,000 

>5.6  X  10* 

R11W12 

>5.6  X  10* 

28,141,500 

>5.6  X  10* 

Notes: 

RCON  =  Concrete  flexural  strength  (Ib/in.^) 

SCI  =  Surface  condition  index  at  failure 

Traffic  simulated  at  rate  of  562,830  ESAL  applications/yr. 


4480  kPa  and  the  SCI  was  increased  to 
90.  A  comparison  of  the  cumulative  dam¬ 
age  with  time  for  the  case  representing 
section  ML5-R-5  with  a  1 .8-m  (6-ft)  wa¬ 
ter  table  (r5w6)  using  all  three  combina¬ 
tions  of  flexural  strength  and  SCI  is  also 
illustrated  in  Figure  14.  Appendix  B  has 
a  compilation  of  graphs  depicting  the 
frost/thaw  penetration  and  cumulative 
damage  for  all  rigid  sections  simulated. 

Cases  with  a  flexural  strength  of 4480 
kPa  exhibited  very  little  damage,  with 
the  MLlO-R-11  section  having  unmea¬ 
surable  damage.  Reducing  the  flexural 
strength  to  3448  kPa  greatly  increased 
the  predicted  damage  for  the  5-year  sec¬ 
tions,  but  had  a  lesser  effect  on  damage 
in  the  10-year  section. 

The  pattern  of  damage  accumulation 


tern  of  damage  accumulation  from  the  Corps  of  with  time  varied  with  the  test  section  simulated 


Engineers  model  is  similar  to  the  MS-1  model, 
especially  when  the  water  table  is  shallow,  caus¬ 
ing  larger  resultant  deflections  and  strains.  For 
simulations  with  deeper  water  tables,  smaller 
amounts  of  damage  accumulate  during  the  sum¬ 
mer. 

Damage  from  models  using  the  vertical  strain 
criteria  follow  three  general  patterns:  1)  a  sharp 
rise  during  spring  thaw  period  with  no  accumula¬ 
tion  at  other  times  of  the  year  (e.g..  Fig.  12),  2)  a 
series  of  increases  during  winter  thaw  events  and 
the  spring  thaw  (Fig.  1 1),  and  3)  in  the  cases  with 
lower  modulus  subgrades,  a  gradual  increase  in 
damage  during  the  spring,  with  sharper  increases 
in  the  summer  (Fig.  13).  This  last  result  is  likely 
due  to  the  anomalously  low  asphalt  moduli  being 
predicted  by  the  Schmidt  model. 

Results — rigid  sections 

Predicted  applications  to  failure  for  the  con¬ 
crete  (rigid)  sections  is  shown  in  Table  14,  with 
the  three  sections  in  the  table  indicating  the  re¬ 
sults  of  a  sensitivity  study.  In  the  left-hand  sec¬ 
tion,  the  flexural  strength  of  the  concrete  was  set 
at  4480  kPa  (650  Ib/in.^)  and  the  surface  condi¬ 
tion  index  (SCI)  at  the  point  of  failure  was  set  at 
80.  In  the  central  section,  the  flexural  strength 
was  reduced  to  3448  kPa  (500  Ib/in.^),  while  the 


(cm)  (in.) 


SCI  remained  at  80.  For  predictions  in  the  right-  Elapsed  ume  (days) 

hand  column,  the  flexural  strength  was  kept  at  Figure  14.  Cumulative  damage  for  case  r5w6. 


(App.  B).  ML5-R-5  experienced  rapid  increases 
in  damage  during  the  thaw  periods,  but  also  had 
an  equivalent  increase  in  damage  over  the  sum¬ 
mer  months.  The  other  sections  experienced  a 
majority  of  their  damage  during  the  thaw  periods 
and  very  little  damage  during  the  summer. 

The  low  damage  predictions  result  from  the 
fact  that  the  model  predicts  damage  at  the  center 
of  a  slab,  whereas  the  majority  of  damage  occurs 
at  the  slab  edges.  Chou  (1989)  recommends  in¬ 
creasing  the  stresses  computed  from  the  layered 
elastic  method  by  a  factor  of  1.33  to  estimate  the 
coverage  levels  for  roadway  pavements. 

Phase  2A 

To  address  problems  discovered  in  the  Phase  1 
modeling  series,  three  additional  simulation  se¬ 
ries  were  conducted  on  the  eight  flexible  test  sec¬ 
tions  studied  in  Phase  1  with  the  water  table  lo¬ 
cated  at  the  shallowest  position.  A  change  that 
applied  to  all  three  series  was  that  the  normalized 
volumetric  unfrozen  water  content  form  of  the 
frozen  modulus  equations  (Table  10a)  was  used, 
rather  than  the  normalized  gravimetric  form  used 
in  Phase  1 . 

The  first  of  the  new  series  was  essentially  a 
repeat  of  the  initial  series  in  that  it  used  the  Schmidt 
model  (eq  8)  for  calculation  of  the  asphalt  modu¬ 
lus  and  the  high  density  1 206  form  of  the  subgrade 
for  modulus  calculations.  It  differed  from  the  ini¬ 
tial  series  in  that  the  normalized  volumetric  unfro¬ 
zen  water  content  form  of  the  frozen  modulus 
equations  was  utilized.  A  mistake  in  the  CUMDAM 
program  was  also  remedied — correcting  the  value 
of  the  subgrade  modulus  used  in  the  Corps  of 
Engineers  vertical  strain  damage  calculation.  Table 
15  lists  the  results  of  this  series  in  the  rows  la¬ 
beled  “Schmidt.”  In  general,  predictions  of  appli¬ 
cations  to  failure  in  the  conventional  cross  sec¬ 
tions  were  less  than  in  the  original  series  using  the 
horizontal  strain  criteria  and  greater  using  the  ver¬ 
tical  strain  criteria.  In  the  full-depth  sections  (F4 
and  F14),  applications  to  failure  were  increased  in 
both  the  horizontal  and  vertical  strain  criteria. 

The  second  series  addressed  the  issue  that  the 
Phase  1  simulations  computed  summer  asphalt 
moduli  that  were  lower  than  considered  reason¬ 
able.  In  this  series,  the  Ullidtz  model  (eq  9)  was 
used  to  calculate  the  asphalt  modulus  at  tempera¬ 
tures  greater  than  1°C;  at  colder  temperatures,  the 


Schmidt  model  was  used.  This  series  also  used  the 
high  density  1206  form  for  the  subgrade  modulus 
calculations.  Results  for  this  series  are  shown  in 
Table  15  in  rows  labeled  “Ullidtz.”  In  nearly  all 
cases,  this  series  produced  higher  predicted  appli¬ 
cations  to  failure  than  the  corresponding  Schmidt 
series. 

The  third  series  also  used  the  Ullidtz  model  to 
calculate  the  asphalt  modulus.  In  addition,  it  ad¬ 
dressed  the  problem  from  the  Phase  1  simulations 
that  during  the  summer  the  predicted  subgrade 
moduli  based  on  the  high-density  1206,  or  “nor¬ 
mal”  form  were  higher  than  FWD-measured  val¬ 
ues.  For  this  series,  the  unfrozen  subgrade  modu¬ 
lus,  when  recovered,  was  set  to  be  constant  at 
103,000  kPa  (15,000  Ib/in.^).  This  value  was  cho¬ 
sen  as  being  approximately  equal  to  an  average 
value  back-calculated  from  FWD  measurements 
on  subgrade  during  fall  1991.*  Asa  result  of  this 
modulus  approximation,  the  linear  model  (model 
0  in  Table  4)  was  used  in  NELAPAV  for  all 
subgrade  calculations.  The  applications  to  failure 
from  this  series  are  listed  in  Table  15  as  “Ull- 
15K.”  In  the  conventional  sections  including  a 
base/subbase,  the  applications  to  failure  were 
slightly  lower  than  those  resulting  from  the  previ¬ 
ous  case  with  the  Ullidtz  asphalt  model  and  the 
higher-modulus  subgrade;  in  the  full-depth  sec¬ 
tions,  they  were  quite  a  bit  lower. 

A  comparison  of  the  Ullidtz/15K  subgrade  se¬ 
ries  results  with  the  originally  designed  failure  at 
3,300,000  applications  for  the  5-year  sections  and 
at  6,600,000  applications  for  the  10-year  sections 
yields  the  following  conclusions: 

1.  In  general,  the  Mechanistic  Pavement  De¬ 
sign  Procedure  predicts  the  5-year  conventional 
sections  to  fail  due  to  asphalt  cracking  prior  to  the 
end  of  the  design  life,  while  failure  from  subgrade 
rutting  is  predicted  beyond  the  design  life. 

2.  Predictions  for  the  5-year  full-depth  section 
indicate  that  it  will  not  fail  from  asphalt  cracking, 
but  two  of  the  three  criteria  for  subgrade  rutting 
indicate  early  failure. 

3.  The  two  10-year  conventional  sections  with 
class  5  and  class  6  special  bases  (F21  and  F22)  are 
predicted  to  fail  from  asphalt  cracking  at  about 
the  end  of  the  design  applications,  while  subgrade 
rutting  criteria  indicate  longer  life. 


*  D.  Van  Deusen,  Mn/ROAD,  pers.  comm.  1992. 
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Table  15.  Predicted  applications  to  failure  (x  1000)  from  simulations  of  flexible  test  sections  run  in  Phase  2A 
series. 


Horizontal  criteria _  _ Vertical  criteria 


Cases 

Asp.  Inst. 
MS-1 

Asp.  Inst. 
MS-11 

Corps  of 
Engineers 

Coetzee/ 

Connor 

Asp.  Inst. 
MS-1 

Corps  of 
Engineers 

FAA 

5-YEAR 

F1W9 

Phase  1 

86 

498 

88 

4,380 

136,942 

5,653 

11,879 

Schmidt 

85 

414 

78 

3,557 

244,709 

3,297 

15,665 

Ullidtz 

651 

935 

218 

17,616 

258,179 

3,430 

17,025 

U11-15K 

617 

863 

201 

15,578 

213,193 

3,133 

18,441 

F2W9 

Phase  1 

522 

1,244 

550 

42,703 

32,477 

1,345 

1,526 

Schmidt 

271 

859 

259 

14,634 

210,011 

2,840 

11,389 

Ullidtz 

970 

1,429 

376 

37,053 

255,832 

3,262 

15,488 

U11-15K 

908 

1,303 

340 

31,925 

117,747 

1,971 

6,649 

F3W8 

Phase  1 

957 

1,706 

958 

92,419 

116,287 

2,653 

9,611 

Schmidt 

533 

1,224 

487 

32,818 

96,046 

1,197 

7,298 

Ullidtz 

1,255 

1,842 

528 

58,812 

87,260 

1,108 

6,381 

U11-15K 

1,225 

1,705 

503 

54,327 

115,098 

1,404 

10,501 

F4W6 

Phase  1 

28,585 

75,751 

90,926 

o 

A 

3,715 

15 

152 

Schmidt 

45,171 

121,039 

175,884 

O 

A 

4,505 

314 

164 

Ullidtz 

65,218 

123,699 

149,292 

>10'^ 

11,952 

317 

246 

U11-15K 

7,055 

26,561 

8,219 

3,752,200 

7,996 

331 

467 

10-YEAR 

F14W50 

Phase  1 

103,272 

485,198 

605,194 

>10^ 

8,602 

10 

251 

Schmidt 

126,479 

721,577 

970,397 

V 

o 

11,687 

417 

336 

Ullidtz 

209,230 

760,581 

852,773 

O 

A 

17,107 

412 

359 

U1M5K 

19,902 

153,779 

40,903 

>10’ 

26,943 

787 

2,232 

F19W18 

Phase  1 

396 

6,148 

966 

75,751 

2,962,263 

56,680 

457,585 

Schmidt 

309 

4,335 

635 

43,328 

4,329,462 

37,175 

632,393 

Ullidtz 

2,948 

14,480 

2,278 

574,316 

4,329,462 

37,422 

639,580 

U11-15K 

2,688 

12,367 

1,976 

469,025 

1,125,660 

18,173 

264,239 

F21W18 

Phase  1 

3,394 

30,456 

12,751 

2,010,107 

562,830 

84,382 

42,478 

Schmidt 

1,958 

16,072 

4,701 

420,022 

865,892 

10,452 

70,442 

Ullidtz 

8,515 

33,402 

9,343 

3,752,200 

792,718 

9,689 

61,849 

U11-15K 

6,661 

23,335 

6,513 

2,345,125 

457,585 

8,752 

66,686 

F22W45 

Phase  1 

7,614 

57,315 

31,147 

8,040,429 

73,190 

22,415 

3,388 

Schmidt 

1,750 

21,507 

5,070 

493,711 

68,388 

1,078 

3,077 

Ullidtz 

8,235 

42,964 

9,192 

3,752,200 

85,667 

1,278 

4,297 

U11-15K 

6,186 

26,045 

5,986 

2,084,556 

120,520 

2,086 

9,840 

Notes: 

Phase  1 — previous  work,  Schmidt —  new  series  with  Schmidt  asphalt  model,  Ullidtz — new  series  with  Ullidtz 
asphalt  model  above  1°C,  U11-15K— new  series  with  Ullidtz  asphalt  model  and  summer  subgrade  modulus  = 
15K  Ib/in .2  (103,000  kPa).  Traffic  simulated  at  rate  of  562,830  ESAL  applications/yr. 
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4.  The  10-year  conventional  section  with  class 
3  special  subbase  (FI 9)  may  fail  due  to  asphalt 
cracking,  but  not  from  subgrade  rutting. 

5.  Predictions  for  the  10-year  full-depth  sec¬ 
tion  indicate  that  it  will  not  fail  from  asphalt  crack¬ 
ing,  but  two  of  the  three  criteria  for  subgrade  rut¬ 
ting  indicate  early  failure. 

Phase  2B 

To  investigate  the  influence  of  the  applied 
freeze/thaw  season  on  the  predicted  performance 
of  the  sections,  an  environmental  effects  sensitiv¬ 
ity  series  was  conducted  using  the  single  conven¬ 
tional  flexible  test  section,  F3,  which  included 
class  5  special  and  class  3  special  materials  as  its 
base  and  subbase  courses  (Fig.  6)  and  a  water 
table  at  a  depth  of  2.4  m  (8  ft).  Prior  simulations 
had  applied  upper  boundary  temperatures  from 
the  season  with  a  freezing  index  near  to  the  30- 
year  mean,  1959-60  (freezing  index  1003  °C-days, 
1806  °F-days).  This  new  series  applied  tempera¬ 
tures  from  three  other  freeze/thaw  seasons,  in¬ 
cluding  those  with  freezing  indices  equal  to  the 
mean  of  the  three  coldest  seasons  during  the  pe¬ 
riod  1959  to  1987  and  the  maximum  and  mini¬ 
mum  for  the  1959-1987  period.  These  seasons 
were  maximum — 1978-79  (1477  °C-days,  2658 
°F-days),  minimum — 1986-87  (467  °C-days,  841 


°F-days),  and  mean  3  coldest  out  of  30 — 1964-65 
(1391  °C-days,  2503  °F-days).  Results  from  these 
years  were  compared  with  F3W8  predictions  for 
the  mean  year  that  had  been  run  in  the  “Schmidt” 
series,  since  this  Phase  2B  series  also  used  the 
Schmidt  model  for  asphalt  modulus  calculations 
(eq  8)  and  the  high  density  1206  form  for  the 
subgrade  modulus  calculations. 

Figure  15  shows  the  freeze/thaw  penetration  in 
the  F3  section  predicted  by  FROST  for  all  four 
freeze  seasons.  In  the  mean  year,  multiple  small 
freeze/thaw  events  were  followed  by  a  long  freeze 
event  that  penetrated  the  subgrade  and  a  spring 
warming  that  remained  above  freezing.  The  mini¬ 
mum  year  had  several  freeze/thaw  events,  but  they 
were  less  severe — none  of  them  penetrated  be¬ 
yond  the  base  course  materials.  The  maximum 
year  began  with  a  very  severe  freeze  event  that 
penetrated  to  the  subgrade  and  lasted  all  winter.  In 
the  spring,  thawing  of  this  freeze  bulb  was  inter¬ 
rupted  by  two  smaller  freeze  events,  such  that 
there  existed  a  thawed  layer  in  the  subbase  be¬ 
tween  two  frozen  zones.  The  year  with  a  freezing 
index  equal  to  the  mean  of  the  three  coldest  in  30 
(3/30  year)  had  a  freeze  season  consisting  of  a 
single  severe  freeze  event  that  penetrated  into  the 
subgrade,  which  then  thawed  with  no  small  freeze 
events. 


(cm)  (in.) 


0  60  120  180  240  300  0  60  120  180  240  300 

Elapsed  Time  (days) 


Figure  15.  Frost  and  thaw  penetrations  predicted  by  FROST  for  freeze 
seasons  in  Phase  2B  environmental  sensitivity  study.  Simulations  start  on 
1  October. 
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Table  16.  Performance  predictions  for  F3  test  section  with  water  table  at  2.4 
m  (8  ft).  Freeze/  thaw  seasons  applied  as  upper  boundary  temperatures  were 
meau  year— 1959-60,  maximum  year— 1978-79,  minimum  year— 1986-87, 
average  of  the  3  coldest  out  of  30 — 1964—65. 


Applications  to  failure  (x  1000) 


Failure  criteria 

Mean 

Max 

Min 

3/30 

Horizontal 

Asphalt  Institute  (MS-1) 

533 

384 

573 

397 

Asphalt  Institute  (MS-11) 

1,224 

1,215 

1,117 

1,445 

Corps  of  Engineers 

A%1 

375 

593 

368 

Coetzee/Conner 

32,818 

25,284 

50,660 

23,530 

Vertical 

Asphalt  Institute  (MS-1) 

96,046 

51,120 

>10^ 

96,046 

Corps  of  Engineers 

1,197 

600 

>10'^ 

1,143 

FAA 

7,298 

4,080 

>10’^ 

7,625 

Note:  Traffic  simulated  at  rate  of  562,830  ESAL  applications/yr. 


Table  16  lists  the  predicted  damage  in  terms  of 
applications  to  failure  from  this  environmental  se¬ 
ries,  in  which  traffic  was  simulated  at  562,830 
ESALs  per  year.  Compared  with  the  mean  year, 
failure  due  to  the  horizontal  strain  criteria  (asphalt 
cracking)  is  predicted  earlier  with  the  maximum 
year  and  the  3/30  year,  and  later  with  the  mini¬ 
mum  year.  However,  the  spread  of  values  is  not 
that  great,  and  3  of  the  4  models  predict  failure 
sooner  than  the  design  figure  of  3,300,000  appli¬ 
cations  for  all  four  of  the  seasons  modeled.  Time 
of  failure  due  to  the  vertical  strain  criteria  (asphalt 
rutting)  is  predicted  sooner  in  the  maximum  year 
than  the  mean  year,  and  at  about  the  same  time 
as  the  mean  year  with  the  3/30  year.  No  damage 
due  to  subgrade  rutting  was  predicted  with  the 
minimum  year  in  which  no  frost  penetrated  the 
subgrade. 

Phase  3 

The  third  phase  of  the  study  modeled  two 
flexible  pavement  sections  to  which  21  years  of 
environmental  conditions  were  applied.  Various 
analyses  were  conducted  attempting  to  correlate 
the  predicted  damage  with  the  characteristics  of 
the  freeze  seasons  simulated. 

Cross  sections/material  properties 

The  two  cross  sections  simulated  employed 
layer  thicknesses  from  two  Mn/ROAD  test  sec¬ 
tions  (Table  7).  The  first,  section  F4,  was  a  “full- 


depth”  section  consisting  of  an  asphalt  layer  lying 
directly  above  the  lean  clay  subgrade.  The  sec¬ 
ond,  section  F3,  was  a  “conventional”  section, 
including  class  6  special  as  the  base  material  (sub¬ 
stituted  for  class  5  special  in  the  actual  pavement 
stmcture)  and  class  3  special  as  the  subbase. 

Physical  properties  used  for  the  materials  were 
the  same  as  shown  in  Table  13,  with  the  exception 
that  a  lower  density  (1.69  Mg/m^;  105.5  Ib/ft^) 
was  used  for  the  1206  subgrade.  Figure  16  illus¬ 
trates  predicted  moduli  for  several  densities  of  the 
1206  subgrade.  The  asphalt  modulus  was  calcu- 


Figure  16.  Resilient  modulus  vs.  degree  of  satura¬ 
tion  of  never-frozen  1206  subgrade  material  illus¬ 
trating  the  effect  of  dry  density. 
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lated  with  a  combination  of  the  Ullidtz  model  at 
above  l^C  temperatures  and  the  Schmidt  model  at 
colder  temperatures. 

Environmental  conditions 

The  series  simulated  the  21 -year  time  period 
1969-70  to  1989-90,  using  a  combination  of  re¬ 
corded  air  temperatures  and  an  estimated  position 
of  the  water  table.  The  simulations  started  on  1 
October  of  the  first  year  and  proceeded  for  365 
days,  using  mean  daily  air  temperatures  recorded 
at  Buffalo,  Minnesota,  applied  as  a  step  condition 
for  24-hr  periods.  The  depth  to  the  water  table 
was  varied  based  on  the  total  precipitation  accu¬ 
mulated  at  St.  Cloud,  Minnesota,  during  the  thaw 
season  prior  to  the  simulated  freeze  season.  An 
inverse  relationship  between  precipitation  and 


water  table  depth  was  applied  such  that  the  maxi¬ 
mum  precipitation  produced  the  shallowest  water 
table  at  0.91  m  (3  ft)  below  the  surface  and  the 
minimum  precipitation  produced  the  deepest  wa¬ 
ter  table  at  1 ,82  m  (6  ft).  The  water  table  was  held 
constant  throughout  the  simulation  at  the  depth 
determined  by  this  relationship.  Figure  17  shows 
the  distribution  of  estimated  water  table  and  freeze 
index  for  this  period. 

Results 

As  in  previous  simulation  series,  there  was  a 
wide  variation  in  the  damage  amounts  predicted 
by  the  different  models  (Table  17).  In  general, 
less  damage  related  to  horizontal  strain  was  pre¬ 
dicted  in  the  full-depth  design  than  in  the  conven¬ 
tional  design.  On  the  other  hand,  when  all  other 


Table  17.  Results  of  21 -year  freezing  index/water  table  series. 


Maximum  Applications  to  failure  (millions) 


Year 

V^ater  table 
depth 
m(ft) 

Freezing 

index 

°C(°F)-days 

Freeze/thaw  freeze  depth 
events  in  into  subgrade 
sub  grade  (cm) 

Horizontal  strain 
Asphalt  cracking 
CEH  AIH-MSll 

Vertical  strain 
Sub  grade  rutting 
CEV  AIV-MSl 

A.  Full  depth  (F4) 

6970  1.5  (4.8) 

1191  (2144) 

2 

91.8 

18.20 

66.00 

0.21 

5.38 

7071 

1.3  (4,2) 

1152  (2074) 

4 

91.8 

18.02 

52.59 

0.34 

8.23 

7172 

1.1  (3.6) 

1274  (2294) 

2 

91.8 

15.16 

45.02 

0.31 

7.64 

7273 

1.2  (4.0) 

839(1510) 

2 

69.8 

18.97 

50.69 

0.37 

7.72 

ITHA 

1.3  (4.2) 

923  (1661) 

2 

77.8 

17.80 

51.48 

0.78 

10.97 

7475 

1.4  (4.6) 

968  (1743) 

3 

81.8 

17.21 

47.93 

0.18 

5.47 

7576 

1.3  (4.2) 

931 (1675) 

3 

73.8 

16.83 

52.26 

0.16 

4.15 

7677 

1.6  (5.3) 

1256  (2261) 

4 

107.8 

17.92 

62.38 

0.16 

4.87 

7778 

1.1  (3.7) 

1331  (2395) 

1 

91.8 

13.88 

60.22 

0.10 

3.01 

7879 

1.2  (3.8) 

1477  (2658) 

3 

99.8 

17.66 

51.36 

6.66 

17.18 

7980 

1.2  (3.8) 

903  (1625) 

5 

79.8 

16.67 

43.10 

0.19 

5.30 

8081 

1.2  (4.0) 

666(1199) 

3 

67.8 

17.15 

42.63 

0.26 

7.72 

8182 

1.8  (6.0) 

1227  (2209) 

3 

107.8 

16.29 

47.92 

0.04 

1.44 

8283 

1.2  (3.9) 

589(1061) 

8 

63.8 

16.73 

34.26 

0.11 

3.71 

8384 

1.0  (3.3) 

1179  (2123) 

2 

81.8 

16.86 

48.96 

0.12 

3.50 

8485 

1.0  (3.3) 

942 (1696) 

3 

81.8 

16.08 

44.68 

0.15 

4.31 

8586 

1.0  (3.2) 

1197  (2154) 

2 

81.8 

17.67 

71.27 

0.79 

11.48 

8687 

0.9  (3.0) 

467  (841) 

4 

45.8 

15.91 

33.52 

0.48 

8.44 

8788 

1.5  (4.8) 

898  (1616) 

4 

87.8 

16.91 

48.28 

0.34 

6.15 

8889 

1.5  (4.9) 

1032(1858) 

4 

87.8 

17.57 

50.53 

0.26 

7.19 

8990 

1.4  (4.7) 

750(1350) 

7 

65.8 

16.50 

37.45 

0.25 

6.83 

Notes: 

CEH  =  Corps  of  Engineers  (horizontal)  CEV  =  Corps  of  Engineers  (vertical) 

AIH-MS 11=  Asphalt  Institute  (MS-1 1)  AIV  =  Asphalt  Institute  (MS-1) 

ND  =  no  damage 

Traffic  simulated  at  rate  of  660,000  ES  AL  applications/yr. 
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Table  17  (cont’d).  Results  of  21-year  freezing  index/water  table  series. 


Maximum  Applications  to  failure  (millions) 


Year 

Water  table 
depth 
m(ft) 

Freezing 

index 

^CCFjdays 

Freeze/thaw  freeze  depth  Horizontal  strain 
events  in  into  subsrade  Asphalt  cracking 
subgrade  (cm)  CEH  AlH-MSll 

Vertical  strain 
Subgrade  rutting 
CEV  AIV-MSl 

A  B.  Conventional  (F3) 

6970 

1.5  (4.8) 

1191  (2144) 

1 

13.4 

0.46 

2.46 

0.82 

61.97 

7071 

1.3  (4.2) 

1152  (2074) 

1 

13.4 

0.56 

2.63 

0.47 

34.83 

7172 

1.1  (3.6) 

1274  (2294) 

1 

17.4 

0.56 

2.53 

0.34 

26.17 

7273 

1.2  (4.0) 

839(1510) 

0 

0 

0.49 

1.69 

ND 

ND 

1314 

1.3  (4.2) 

923  (1661) 

2 

1.4 

0.46 

2.04 

ND 

ND 

7475 

1.4  (4.6) 

968  (1743) 

1 

5.4 

0.49 

2.61 

1.46 

120.66 

7576 

1.3  (4.2) 

931  (1675) 

1 

0 

0.55 

2.71 

ND 

ND 

7677 

1.6  (5.3) 

1256  (2261) 

1 

25.4 

0.43 

2.18 

0.71 

58.36 

111% 

1.1  (3.7) 

1331  (2395) 

1 

13.4 

0.49 

2.31 

0.51 

38.69 

im 

1.2  (3.8) 

1477  (2658) 

1 

21.4 

0.49 

2.72 

0.36 

27.68 

7980 

1.2  (3.8) 

903  (1625) 

1 

5.4 

0.52 

1.87 

0.44 

33.62 

8081 

1.2  (4.0) 

666(1199) 

0 

0 

0.50 

1.50 

ND 

ND 

8182 

1.8  (6.0) 

1227  (2209) 

1 

29.4 

0.42 

2.21 

0.46 

37.95 

8283 

1.2  (3.9) 

589(1061) 

0 

0 

0.51 

1.34 

ND 

ND 

8384 

1.0  (3.3) 

1179(2123) 

1 

5.4 

0.49 

1.84 

1.18 

91.16 

8485 

1.0  (3.3) 

942  (1696) 

1 

5.4  . 

0.48 

1.56 

0.98 

77.74 

8586 

1.0  (3.2) 

1197  (2154) 

1 

5.4 

0.48 

2.70 

1.64 

ND 

8687 

0.9  (3.0) 

467  (841) 

0 

0 

0.44 

1.45 

ND 

ND 

8788 

1.5  (4.8) 

898(1616) 

1 

9.4 

0.45 

1.85 

0.72 

58.30 

8889 

1.5  (4.9) 

1032  (1858) 

1 

9.4 

0.53 

1.89 

0.45 

33.69 

8990 

1.4  (4.7) 

750  (1350) 

0 

0 

0.55 

1.65 

ND 

ND 

Notes:  CEH  =  Corps  of  Engineers  (horizontal)  CEV  =  Corps  of  Engineers  (vertical) 

AIH-MS 1 1  =  Asphalt  Institute  (MS-1 1)  AIV  =  Asphalt  Institute  (MS-1) 


ND  =  no  damage 

Traffic  simulated  at  rate  of  660,000  ESAL  applications/yr. 


(m)  (ft) 


Freezing  Index 


Figure  17.  Distribution  of  freezing  indices  and  water  table  depths  in  21 -year 
Phase  3  series. 
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conditions  are  constant,  there  was  more  damage 
related  to  vertical  strain  in  the  full-depth  design. 
When  frost  did  not  penetrate  the  subgrade  in  the 
conventional  design,  essentially  no  vertical  strain 
damage  was  predicted. 

Freeze  season  characteristics  were  initially 
quantified  by  analyzing  the  total  freezing  index, 
the  number  of  times  that  frost  was  predicted  to 
enter  the  subgrade  and  the  maximum  depth  of 
frost  penetration  beneath  the  top  of  the  subgrade 
(Table  17).  Lower  amounts  of  damage  (higher 
applications  to  failure)  were  predicted  by  the  As¬ 
phalt  Institute  horizontal  strain  model  in  seasons 
with  higher  freezing  indices  and  deeper  frost  pen¬ 
etration  in  both  the  full-depth  and  conventional 
cross  sections.  However,  predictions  of  damage 


Figure  18.  Predicted  results  from  simulation  of 
fulUdepth  section  during  freeze  season  1983- 
1984,  starting  on  1  October,  a)  Frost/thaw  depth, 
b)  seed  resilient  modulus  passed  to  NELAPAV, 
and  c)  cumulative  damage.  (Abbreviations — see 
Table  17.). 


from  the  other  models  did  not  correlate  with  any 
of  the  above  parameters. 

It  was  noted  that  two  similarly  categorized  sea¬ 
sons,  1983-84  and  1985-86,  which  had  nearly 
identical  freezing  indexes  and  water  table  depths, 
as  well  as  the  same  number  of  freeze  thaw  events, 
had  diversely  different  damage  predictions  in  the 
full-depth  section  (Table  17a).  A  large  thaw  event 
in  the  middle  of  the  1983-1984  freeze  season  was 
accompanied  by  severe  weakening  of  the  subgrade 
both  during  midwinter  and  spring  thaw,  which 
resulted  in  a  considerable  amount  of  predicted 
damage  (Fig.  18).  The  1985-86  year  experienced 
a  nearly  continuous  freeze  season  with  no  signifi¬ 
cant  midwinter  thaw  event  (Fig.  19).  As  a  result, 
the  subgrade  modulus  underwent  thaw  weakening 


Figure  19.  Predicted  results  from  simulation  of 
full-depth  section  during  freeze  season  1985- 
1986,  starting  on  1  October,  a)  Frost/thaw  depth, 
b)  seed  resilient  modulus  passed  to  NELAPAV, 
and  c)  cumulative  damage  (see  Table  17). 
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Annual  Damage 


only  during  spring  thaw  in  the  1985-86  season,  and 
the  corresponding  total  damage  accumulation  was 
much  less. 

Based  on  the  above  observation,  we  tried  a  new 
approach  to  quantifying  the  characteristics  of  the 
freeze  seasons.  It  involved  summing  the  thawing 
degree  days  experienced  during  the  freeze  season 
in  two  ways.  The  first  was  a  total  from  all  events, 
and  the  second  was  a  total  from  only  the  thaw 
events  that  exceeded  16°C-days  (30°F-days)  (Table 
18).  The  quantity  16°C-days  was  chosen  based  on 
analysis  by  Mahoney  et  al.  (1985)  indicating  that 
pavements  approach  their  critically  weak  condi¬ 
tion  after  this  amount  of  thaw  has  been  experi¬ 
enced.  Unfortunately,  neither  of  these  quantities 
correlated  with  predicted  damage  amounts.  Ap¬ 


parently,  a  more  sophisticated  index  that  perhaps 
combines  freeze  index,  severity  of  midwinter  thaw 
events,  and  other  parameters  is  required  to  corre¬ 
late  with  predicted  damage. 

In  analyzing  the  Phase  3  simulation  results,  we 
examined  the  distribution  of  predicted  damage 
through  the  four  seasons  defined  as  follows:  fall — 
1  October  to  start  of  freeze,  winter — freeze  sea¬ 
son  as  defined  by  freezing  index,  spring — 75  days 
following  end  of  freeze  season,  summer — remain¬ 
der  of  365-day  simulation.  The  simulations  pre¬ 
dicted  a  wide  variation  in  seasonal  damage 
amounts  for  different  years,  especially  in  the  ver¬ 
tical  strain  damage  in  the  conventional  cross  sec¬ 
tion  (Fig.  20).  In  years  when  vertical  strain  dam¬ 
age  was  predicted  for  the  conventional  cross 


a.  Conventional  section.  Asphalt  Institute  horizontal  c.  Full  depth  section,  Asphalt  Institute  horizontal  strain 
strain  model  MS-1 1.  model  MS-1 1 . 


b.  Conventional  section.  Asphalt  Institute  vertical  strain  d.  F ull  depth  section.  Asphalt  Institute  vertical  strain 
model  MS-1.  model  MS-1. 


Figure  20.  Distribution  of  cumulative  damage  during  seasons. 


Table  18.  Midwinter  thaw  index  summations. 


Water  table 

Freezing 

Midwinter  thaw  index 

depth 

index 

(Total) 

(>  30  DD) 

Year 

m(ft) 

^CCF)-days 

°C(°F)-days 

°C(°F)-days 

6970 

1.5  (4.8) 

1191  (2144) 

25.5  (45.9) 

0.0 

7071 

1.3  (4.2) 

1152  (2074) 

41.0  (73.8) 

0.0 

7172 

1.1  (3.6) 

1274  (2294) 

50.0  (90.0) 

19.4 

(35.0) 

7273 

1.2  (4.0) 

839  (1510) 

21.5  (38.7) 

0.0 

7374 

11.3  (4.2) 

923  (1661) 

56.0  (100.8) 

22.3 

(40.1) 

7475 

1.4  (4.6) 

968  (1743) 

35.0  (63.0) 

14.2 

(25.5) 

7576 

11.3  (4.2) 

931  (1675) 

67.0  (120.6) 

20.7 

(37.2) 

7677 

1.6  (5.3) 

1256  (2261) 

44.0  (79.2) 

0.0 

7778 

1.1  (3.7) 

1331  (2395) 

26.5  (47.7) 

0.0 

7879 

1.2  (3.8) 

1477  (2658) 

52.0  (93.6) 

18.0 

(32.4) 

7980 

1.2  (3.8) 

903  (1625) 

32.0  (57.6) 

0.0 

8081 

1.2  (4.0) 

666  (1199) 

41.0  (73.8) 

17.5 

(31.5) 

8182 

1.8  (6.0) 

1227  (2209) 

31.5  (56.7) 

21.9 

(39.4) 

8283 

11.2  (3.9) 

589  (1061) 

92.5  (166.5) 

75.0  (135.0) 

8384 

1.0  (3.3) 

1179  (2123) 

40.0  (72.0) 

34.7 

(62.4) 

8485 

1.0  (3.3) 

942  (1696) 

40.0  (72.0) 

22.9 

(41.2) 

8586 

1.0  (3.2) 

1197  (2154) 

29.0  (52.2) 

0.0 

8687 

10.9  (3.0) 

467  (841) 

65.0  (117.0) 

40.1 

(72.1) 

8788 

11.5  (4.8) 

898  (1616) 

79.0  (142.2) 

54.2 

(97.6) 

8889 

1.5  (4.9) 

1032  (1858) 

49.0  (88.2) 

0.0 

8990 

11.4  (4.7) 

750  (1350) 

70.5  (126.9) 

27.0 

(48.6) 

Table  19.  Average  percentage  of  total  yearly  dam¬ 
age  accumulated  during  four  seasons. 


Accumulated  damage  (%  of  total) 


Section/model 

Fall 

Winter 

Spring 

Summer 

Conventional 

AIH 

22 

31 

45 

2 

AIV 

0 

40 

60 

0 

CEV 

0 

42 

58 

0 

Full  depth 

AIH 

16 

51 

32 

1 

AIV 

2 

44 

16 

38 

CEV 

0 

80 

20 

0 

Notes: 

AIH  =  Asphalt  Institute  horizontal  (MS-1 1) 

AIV  =  Asphalt  Institute  vertical  (MS-1) 

CEV  =  Corps  of  Engineers  vertical 

section,  a  majority  occurred  in  the  spring,  while 
the  remainder  occurred  in  the  winter  (Table  19). 
Horizontal  strain  damage  in  the  conventional  cross 
section  also  occurred  mainly  in  the  spring,  with 
additional  significant  amounts  in  both  the  winter 
and  fall,  and  very  small  amounts  in  the  summer. 


In  the  full-depth  cross  section,  horizontal  strain 
damage  occurred  primarily  in  the  winter,  with 
slightly  less  in  the  spring,  and  some  additional 
damage  in  the  fall.  Both  vertical  strain  models 
predicted  damage  in  the  full-depth  section  to  oc¬ 
cur  mainly  in  the  winter,  with  some  in  the  spring. 
The  Asphalt  Institute  vertical  strain  model  also 
predicted  additional  amounts  of  summer  damage. 

DISCUSSION  AND 
RECOMMENDATIONS 

The  Phase  1  modeling  series  indicated  signifi¬ 
cantly  different  performance  by  the  different  test 
sections  and  highly  variable  results  depending  on 
the  performance  model  applied  (Table  13).  The 
simulated  performance  of  the  test  sections  was 
also  significantly  affected  by  the  subgrade  condi¬ 
tions,  e.g.,  density,  soil  moisture  and  water  table 
depth.  For  example,  compare  the  model  predic¬ 
tions  using  the  Asphalt  Institute  MS-1  horizontal 
strain  criteria.  For  case  f4w6,  the  full-depth  5- 
year  section  with  the  1206  subgrade  in  its  opti- 
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mum  density  condition,  the  model  predicts 
28,585,000  applications  to  failure.  In  case  f4w61d, 
the  same  test  section  with  a  low  density  value  for 
the  1206  subgrade,  the  model  predicts  1,905,000 
applications  to  failure.  And  finally,  in  case  f4w6ss, 
with  the  1232  subgrade  that  produced  lower  resil¬ 
ient  modulus  values  in  lab  tests,  failure  is  pre¬ 
dicted  after  only  72,000  applications,  or  in  less 
than  one  year. 

After  the  Phase  1  simulations  were  complete, 
the  resilient  modulus  data  from  the  1206  subgrade 
were  reviewed.  A  calibration  error  of  unknown 
magnitude  was  discovered  in  the  resilient  modu¬ 
lus  equipment  used  to  measure  the  unfrozen  modu¬ 
lus  on  that  material.  Modulus  testing  of  the  1206 
subgrade  in  the  frozen  condition  and  all  testing  on 
the  other  materials  were  conducted  on  a  different 
machine,  which  passed  its  calibration  checks.  Com¬ 
paring  frozen  and  unfrozen  resilient  modulus  test 
results  from  the  1206  subgrade  with  those  from 
other  materials  tested.  Figure  7,  indicates  that  the 
unfrozen  values  for  the  1206  subgrade  may  be 
about  an  order  of  magnitude  too  high.  A  lower 
for  the  1206  subgrade  would  have  resulted  in  ear¬ 
lier  failures  than  computed  in  this  report  where 
the  1206  optimum  density  subgrade  was  used. 

We  are  currently  investigating  the  unfrozen 
1206  subgrade  results  and  will  include  findings 
and  revised  performance  predictions  in  Berg  (in 
prep.) 

It  is  obviously  extremely  important  to  use  the 
representative  subgrade  conditions  in  the  design 
simulations.  Results  are  consistent  with  observa¬ 
tions  from  in-service  pavements;  e.g.,  weak  areas 
fail  much  more  rapidly  than  strong  ones  and  high 
water  tables  cause  failures  before  similar  pave¬ 
ments  constructed  over  lower  water  tables.  Both 
the  1206  and  1232  subgrade  samples  were  ob¬ 
tained  from  the  test  site  and  were  located  less  than 
one-half  mile  apart,  yet  their  tested  behavior  was 
quite  different. 

Another  important  aspect  governing  the  results 
of  the  Mechanistic  Pavement  Design  Procedure  is 
the  method  used  to  calculate  the  asphalt  concrete 
resilient  modulus.  The  use  of  the  Schmidt  model, 
which  produces  extremely  low  summer  season 
moduli,  partially  accounts  for  the  short  lifespans 
predicted  in  the  Phase  1  modeling  series.  Replace¬ 
ment  with  the  Ullidtz  model  in  part  of  the  Phase 
2A  modeling  series  increased  the  predicted 


lifespans  of  the  sections  when  judged  by  the  hori¬ 
zontal  strain  criteria  (Table  15). 

Results  from  the  Phase  2A  series  that  used  the 
Ullidtz  model  for  asphalt  modulus  calculations 
and  a  recovered  summer  subgrade  modulus  of 
103,00  kPa  (15,000  Ib/in.^),  a  close  approxima¬ 
tion  of  the  value  backcalculated  from  FWD  mea¬ 
surements  on  subgrade  during  fall  1991,  yield  per¬ 
formance  predictions  that  are  as  “fine-tuned”  as 
possible,  so  far  (Table  15  “U11-15K”).  When  these 
Mechanistic  Pavement  Design  Procedure  predic¬ 
tions  are  compared  with  the  originally  designed 
failure  at  3,300,000  applications  for  the  5-year  sec¬ 
tions  and  6,600,000  applications  for  the  10-year 
sections,  the  following  statements  can  be  made: 

1.  Predictions  for  the  full-depth  sections,  both 
5  and  10  year,  indicate  that  they  will  not  fail  from 
asphalt  cracking,  but  two  of  the  three  criteria  for 
subgrade  rutting  indicate  early  failure; 

2.  Conventional  sections  are  predicted  not  to 
fail  due  to  subgrade  rutting;  however,  sections 
with  more  frost-susceptible  bases  are  predicted  to 
fail  because  of  asphalt  cracking  relatively  early 
in  their  design  life,  and  sections  with  non-frost- 
susceptible  bases  are  predicted  to  fail  towards  the 
end  of  their  design  life. 

By  modeling  21  years  of  environmental  condi¬ 
tions  in  the  Phase  3  series,  we  were  hoping  to  find 
the  characteristics  of  a  “most  severe”  winter  that 
could  be  used  for  design  purposes.  The  results 
were  complex  and  a  relation  between  predicted 
failure  and  characteristics  of  freeze  seasons  was 
elusive,  at  best. 

Based  on  the  above  results,  we  recommend  the 
following  studies: 

1 .  The  class  4  special  and  class  5  special  base 
materials  should  receive  the  full  complement  of 
laboratory  tests  so  that  simulations  may  be  run 
using  properties  of  actual  materials  from  the  site 
rather  those  of  substitute  materials. 

2.  Once  performance  data  are  received  from 
Mn/ROAD,  predicted  and  measured  values  of 
moisture,  temperature,  and  strains  should  be  com¬ 
pared. 

3.  Use  performance  data  from  Mn/ROAD  to 
revise,  refine  or  develop  new  damage  models  for 
rutting  and  fatigue  cracking. 

4.  Use  performance  data  from  Mn/ROAD  to 
refine  or  develop  a  new  model  for  the  change  of 
the  asphalt  modulus  with  temperature. 
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5.  Use  data  from  falling  weight  deflectometer 
tests  on  Mn/ROAD  to  estimate  pavement  perfor¬ 
mance  on  an  annual  basis. 

6.  The  Procedure  should  be  used  to  estimate 
performance  of  all  of  the  test  sections.  This  study 
evaluated  less  than  one-half  of  the  sections. 


CONCLUSIONS* 

The  range  of  values  produced  for  the  various 
scenarios  is,  as  noted,  extremely  wide.  This  leads 
to  the  conclusion  that  mechanistic  design  in  its 
present  stage,  while  a  powerful  predictor  of 
changes  in  pavement  response  with  changes  in 
loads  or  moduli,  is  at  present  an  uncertain  predic¬ 
tor  of  pavement  performance.  Data  and  analyses 
from  Mn/ROAD  are  cmcial  to  improving  that  cur¬ 
rent  uncertain  state  of  performance  prediction. 
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